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Abstract 
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ysis is performed in the fixed-flavor number scheme for = 3,4,5 and uses the MS -scheme for 
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results are used to compute benchmark cross sections at hadron colliders to NNLO accuracy and 
to compare to data from the LHC. 
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1 Introduction 



Parton distribution functions (PDFs) in the nucleon are an indispensable ingredient of modern 
collider phenomenology and their study has a long history. In the perturbative approach to the 
gauge theory of the strong interactions, Quantum Chromodynamics (QCD), factorization allows 
for the computation of the hard parton scattering processes as a power series in the strong coupling 
constant as and, typically, to leading power l/Q^ dominating for large momentum transfer Q^. 
Predictions for physical cross sections involving initial hadrons, however, do require further non- 
perturbative information, that is knowledge of the PDFs in the nucleon as well as the value of as(Q) 
and of the masses of heavy quarks. Since PDFs cannot be calculated in perturbative QCD, they 
need to be extracted from a comparison of theory predictions to available experimental precision 
data on deep-inelastic scattering (DIS), on the production of lepton-pairs (Drell-Yan process) or 
jets in hadron collisions or any other suitable hard scattering reaction. 

The accuracy of PDF determinations in such analyses has steadily improved over the years, 
both due to more accurate experimental input and due to refined theory predictions for the hard 
parton scattering reactions including higher orders in perturbation theory. As of now, complete 
next-to-next-to-leading order (NNLO) calculations in perturbative QCD form the backbone of this 
endeavor. These allow for the computation of many important benchmark cross sections, e.g. in the 
proton-proton collisions at the LHC, with an unprecedented precision. The determination of PDFs 
to NNLO accuracy in QCD was pioneered more than a decade ago in [1] and builds in particular 
on the known corrections for the PDF evolution |l2l|3l as well as on the hard scattering corrections 
for DIS [|4]-[T0l. and hadronic W- and Z-gauge-boson production, both at an inclusive [1II1II21 and 
a differential level [fBUTSl. 

Presently, NNLO PDFs have been obtained by a number of different groups. In detail, these 
are ABKM09 G61II3, HERAPDF1.5 JR09 Mi^, MSTW [22] and NN21 [|23l, while 

CTIO [|24l still remains at next-to-leading order (NLO) accuracy only. There exist, of course, 
differences between these PDF sets. These arise from variations in the choice of the parameters, 
e.g., the value of OsiMz), but also from a different theoretical footing for the data analysis. In the 
latter case, this comprises for instance, the treatment of the heavy-quark contributions in DIS, the 
corrections for nuclear effects, the inclusion of higher twist (HT) terms and so on. The implications 
for precision predictions at TeV-scale hadron colliders can be profound, though, as benchmark 
cross sections at NNLO in QCD for the production of W- and Z-gauge-bosons or the Higgs boson 
through gluon-gluon-fusion (ggF) show, see e.g., the recent discussion in [23'-'30'|. 

In this article we present the PDF set ABMl 1, which is an updated version of the PDF analyses 
of ABKM09 [m and ABMIO [[13 in the 3-, 4-, and 5-flavor scheme at NNLO in QCD. These 
PDFs are obtained from an analysis of the world DIS data combined with fixed-target data for the 
Drell-Yan (DY) process and for di-muon production in neutrino-nucleon DIS. In the ABMll fit 
we are now using the final version of the DIS inclusive data collected by the HERA experiments 
in run I [[TSlI together with new data of the HI collaboration from the HERA low-energy run OTI . 
Moreover, our update is based on theoretical improvements. For instance, the treatment of the 
heavy-quark contributions in DIS now employs the running-mass definition in the MS -scheme for 
the heavy quarks ll32ll . 

The strong coupling constant as(Mz) or, respectively, the QCD scale Aqcd. is a mandatory 
parameter to be fitted in DIS analyses of world data, its correct value being of paramount impor- 
tance for many processes in DIS and at hadron colliders, in particular for Higgs boson production 
in ggF [26J. An essential criterion for the selection of additional precision data on top of the world 
DIS data, e.g. those for the DY process or for hadronic weak-boson and jet-production cross sec- 
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tions, in the measurement of the QCD scale is the compatibility of these data sets with respect 
to the experimental systematics in the different measurements. Strictly speaking, combined anal- 
yses require a theoretical description at the same perturbative order. Because of these reasons, 
the combination of different data sets needs great care if performed with the goal of a precision 
measurement of a^(Mz). Combinations of a wide range of hard-scattering data sets of differing 
quality, as sometimes used in more global fits, are useful only, if they indeed lead to a statistically 
and systematically improved value of UsiMz). Of course, a careful check is always required when 
new data sets are added. As a result of our new analysis we determine in the ABMl 1 fit the strong 
coupling constant at NNLO in the MS -scheme and present a detailed discussion of the uncertain- 
ties and of the impact of individual experiments, showing the great stability in the obtained value 
of a^iMz). 

The papers is organized as follows. In Sec.[2]we describe the theoretical framework of our anal- 
ysis, in Sec. I2.1l in particular the perturbative QCD input including the framework for heavy-quark 
DIS. Sees. 12.21 and [23] are concerned with a detailed account of the non-perturbative corrections 
and nuclear corrections, which have already been applied in previous PDF analyses ni6[|17L Sec. [3] 
features in detail the data analysed with an emphasis on the systematic and normalization uncer- 
tainties. This comprises data on inclusive DIS from the HERA collider and fixed target experiments 
in Sec. 13. 1[ on the DY process in Sec. 13.21 and on di-muon production in neutrino-nucleon DIS in 
Sec. [331 

The main results of the present work are contained in Sec. [H where we present all PDF param- 
eters along with illustrations of the shapes of PDFs. The numerous checks include studies of the 
pulls and the statistical quality for all individual experiments as well as a detailed assessment of the 
power corrections induced by the higher twist terms. As we work in a scheme with a fixed number 
Hf of light quark flavors, a detailed discussion is also devoted to the generation of heavy-quark 
PDFs. Our determination of the strong coupling constant as at NNLO in QCD leads to the value 
(XsiMz) = 0.1 134 + 0.001 1. We show the impact of the individual data sets on as(Mz) and compare 
with the determinations from other PDF fits and other measurements included in the current world 
average. Finally, Sec. [4] is complemented with a comparison of moments of PDFs with recent 
lattice results. 

The consequences of the new PDF set ABMll on standard candle cross section benchmarks 
are illustrated in Sec.[51 We provide cross section values for W- and Z-boson production in schemes 
with nf = 4 and = 5 flavors and we address the accuracy of theory predictions for all dominant 
Higgs boson search channels at the LHC. The PDF uncertainties for top-quark pair-production 
are also illustrated highlighting the combined uncertainty in the gluon PDF, as and the top-quark 
mass mt. Sec. [5] finishes with comments on the issue of hadronic jet production, especially from 
the Tevatron, and the impact of its data on PDF fits. We conclude in Sec. [6] and summarize our 
approach for the handling of the correlated systematic and normalization uncertainties along with 
the explicit tables for the covariance matrix of the ABMl 1 fit in App.[A[ 

2 Theoretical framework 

Here, we briefly recall the theoretical basis of our PDF analysis, which is conducted in the so- 
called fixed-flavor number scheme (FFNS) for nf light (massless) quarks. That is to say, we con- 
sider QCD with n f light quarks in the PDF evolution, while heavy (massive) quarks only appear 
in the final state. As far as QCD perturbation theory is concerned, we specifically focus on aspects 
relevant to NNLO accuracy. For completeness our treatment of power corrections and also of nu- 
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clear corrections as needed e.g., for DIS from fixed-target experiments, is documented in Sec. 12.21 
and l2.3[ The latter have already been used in our previous PDF determinations [[T6l[T7l . 



2.1 Perturbative QCD 

The ability to make quantitative predictions in QCD which is a strongly coupled gauge theory, 
rests entirely on its factorization property. A cross section for the production of some final state X 
from scattering of initial state hadrons can be expressed in lepton-nucleon {ep) DIS as, 

cTep-^ix = ^ J dzfi(z,as(p,),iuj)d-ei-,x(z,Q^,as(ju^),p.^,p.j) , (2.1) 
for l = e,v and in proton-proton collisions (pp) as, 

o-pp^x = ^ I dzidz2fi(zuas(Mr)'l^f)fj(z2,as(Mr)'l^f)^ 

ij ^ 

xd-ij^x(zi,Z2,Q^,as(Mr)'l^r^l^f) ' (2.2) 
where the PDFs in the nucleon fj (i = q,q,g) are the objects of our primary interest. They describe 



the nucleon momentum fraction z (or zi, zi) carried by the parton and the sums in eqs. (|2.1I) 
and (12.21 ) run over all light (anti-)quarks and the gluon. The parton cross sections denoted & are 
calculable in perturbation theory in powers of the strong coupling constant Us and describe the hard 
interactions at short distances of order (9(1/2). We have also displayed all implicit and explicit 
dependence on the renormalization and factorization scales, p^ and /i^. Throughout our analysis, 
however, we will identify them, p^ = = p. AH dependencies of (jep^ix and cTpp^x on the 
kinematics and, likewise the integration boundaries of the convolutions, has been suppressed in 
eqs. (12.11 ) and (12.21 ), as these are specific to the observable under consideration. 

In standard DIS, the (semi-)inclusive cross section (Tep->ix in eq. (12.11 ) depends on the Bjorken 
variable x, the inelasticity y and on Q^, the (space-like) momentum-transfer between the scattered 
lepton and the nucleon. Moreover, it admits a decomposition in terms of the well-known DIS 
(unpolarized) structure functions F,, / = 1,2,3. QCD factorization applied to the DIS structure 
functions implies 

1 

akFk{x,Q^) = 2 J^fi(z,as(p),l^^)Ck,i[^,Q^,as{p),p^^, (2.3) 

where ai = 2, a2 = l/x, aj, = 1, and Ck,i denote the Wilson coefficients. Fi = F2- Fj defines 
the longitudinal structure function in terms of the tranverse structure function Fj = 2xF\, see also 
eq. (13.21) below for the relation including target masses. Eq. (12.31) integrated over x gives rise to the 
standard Mellin moments, 

1 

Fk{N,Q^) = J dxx^-^akF,(x,Q^). (2.4) 


These link the theoretical description of DIS to the operator-product expansion (OPE) on the light- 
cone. The OPE allows to express the DIS structure functions as a product of (Mellin moments 
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of) the Wilson coefficients Ckj and operator matrix elements (OMEs) of leading twist (twist-2). 
Moreover, it admits a well-defined extension in powers of l/Q^ (twist-4, twist-6 and so on), cf. 
Sec.O 

The scale dependence of the PDFs is contained in the well-known evolution equations 



(2.5) 



at leading twist, which is a system of coupled integro-differential equations corresponding to the 
different possible parton splittings. The splitting functions Pij in eq. (12.51 ) have been determined 
at NNLO in |!2l.'3l, which implies knowledge on the first three terms in the powers series in Os 
(suppressing parton indices). 



a[P^'^ . (2.6) 



The PDFs fi are subject to sum rule constraints due to conservation of the quark number and the 
momentum in the nucleon, which imply at each order in perturbation theory a vanishing first (sec- 
ond) Mellin moment for specific (combinations of) splitting functions Pij in eq. (12.61 ). These sum 
rule constraints relate the PDF fit parameters used in the parametrizations of the input distributions, 
see Sec. |4l The accuracy of the numerical solution of the differential eq. (12.51) up to NNLO was 
tested by comparison to programs such as QCD- PEGASUS [|33]| or HOPPET ||34]| . 

For the massless DIS structure functions we will be using the following input from perturbative 
QCD at leading twist, 

oo oo 

Fk = J] ^'sFf , k = 2,3, FL = a,,J] (^[pf > (2-7) 

/=o /=o 

where, again, (NLO) NNLO accuracy is defined by the first (two) three terms in the power series 
in as , cf. eq. (12.61) . The Wilson coefficients for F2, Fj, in eq. (12.71) are known to NNLO from 
[|4]-[8l, and, actually, even to next-to-next-to-next-to-leading order (N^LO) from [[I0l|35l, and for 
Fl to NNLO from [|9l[l0|. Note that in the latter case the perturbative expansion starts at order 
as, thus NNLO accuracy for Fl actually requires three-loop information, which is numerically not 
unimportant. 

Likewise, for the partonic cross sections of the DY process in eq. (12.21) . i.e., for hadronic W- 
and Z-boson production, we use 

00 

/=0 

with the NNLO resuffs of [fTTHTSH. 

At the level of NNLO accuracy, QCD perturbation theory is expected to provide precise predic- 
tions as generally indicated by the numerical size of the radiative corrections at successive higher 
orders and their pattern of apparent convergence. The residual theoretical uncertainty from the 
truncation of the perturbative expansion is conventionally estimated by studying the scale stability 
of the prediction, i.e., by variation of the renormalization and factorization scales fi^ and //^ in 
eqs. (12.11) and (12.21 ). As stated above, we set //^ = /i^ = /i in our analysis, and moreover, identify the 
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scale jj. with the relevant kinematics of the process, e.g., fi = Q for DIS. Currently no PDF fits with 
an independent variation of /i,. and ju^ are available and we leave this issue for future studies. 

One important aspect is the production of heavy quarks in DIS both for the neutral-current 
(NC) and the charged-current (CC) exchange. In the former case, pair-production of charm-quarks 
accounts for a considerable part of the inclusive DIS cross section measured at HERA, especially 
at small Bjorken-x, while the latter case is needed in the description of neutrino-nucleon DIS. 
At not too large values of Q^, the NC reaction is dominated by the photon-gluon fusion process 
'y*g ^ ccX, while the CC case proceeds through W*s ^ c, so that the perturbative expansion of 
the respective heavy-quark structure functions reads, 

oo 00 

Fl^^{x,Q\ml) = a, J]a[Ff\ Fl^^{x,Q\ml) = J]a[Ff\ (2.9) 

/=0 ' /=0 

where k = 2,3,L and niq is the heavy-quark mass. The heavy-quark Wilson coefficients are known 
exactly to NLO, both for NC [|36l and CC [|371|38l- The NNLO results for F^^^ are, at present, ap- 
proximate only and based on the logarithmically enhanced terms near threshold Il39] - |4TI (see [|42l 
for threshold resummation in the CC case). As well known [|43l . the heavy-flavor corrections to 
F2 are represented with an accuracy of 0(1%) and better for Q^/rn^ > 10. Under this condition the 
Wilson coefficients are given by Mellin convolutions of massive OMEs [|43]-l46l and the massless 
Wilson coefficients [|4]-[T0l. Fixed Mellin moments of the heavy-quark OMEs have also been com- 
puted at three loops in BTl and first results for general values of Mellin-A'^ have been calculated 
in [|48]| . Mellin space expressions for the NC and CC Wilson coefficients up to 0{a^) are available 

in m\M. 

In the current PDF analysis, the bulk of data from DIS experiments can be described in a 
scheme with n/ = 3 light flavors. At asymptotically large scales Q » nicinh the genuine contribu- 
tions for heavy charm- and bottom-quarks in a FENS with = 3 grow as a siQ)ln{Q^ / m^), as the 
quark masses screen the coUinear divergence. The standard PDF evolution equations in eq. (|2.5I) 
resum these logarithms at the expense of matching the eff"ective theories, i.e. QCD with nf and 
nf+l light flavors. This defines a variable-flavor number scheme (VENS) and gives rise to the so- 
called heavy-quark PDEs for charm- and bottom-quarks in QCD with effectively nf = 4 and nf = 5 
light flavors. The heavy-quark PDEs are generated from the light flavor PDEs in a nf = 3-flavor 
FENS as convolutions with OMEs, see e.g. [fT6ll44l . The VENS requires the matching conditions 
both for the strong coupling as and the PDEs (through the corresponding OMEs), which are known 
to N^LO [I5n.(52| for or,, and to NNLO for the OMEs [.44. 46.1 . An extensive discussion of the VENS 
implementation has been presented in our previous analysis [fT6l . 

The heavy-quark masses in eq. (12.91 ) are well-defined within a specific renormalization scheme, 
the most popular ones being the on-shell and the MS -scheme. The former uses the so-called 
pole-mass niq, defined to coincide with the pole of the heavy-quark propagator at each order in 
perturbative QCD, and known to have intrinsic theoretical limitations. As a novelty of our analysis, 
we employ the MS -scheme for niq, which enters both in the massive OMEs and in the Wilson 
coefficients and introduces a running mass mq(ju) depending on the scale of the hard scattering 
in complete analogy to the running coupling a sip). As a benefit, predictions for the heavy-quark 
structure functions in terms of the MS -mass display better convergence properties and greater 
perturbative stability at higher orders ^2\, thus reducing the inherent theoretical uncertainty. 

The Fortran code OPENQCDRAD for the numerical computation of all hard scattering cross 
sections within the present PDF analysis is publicly available [53]. It comprises in particular the 
theory predictions for the DIS structure functions including the heavy-quark contributions as well 
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as for the hadronic W- and Z-boson production and it is capable of computing of the benchmark 
cross sections to NNLO accuracy in QCD in Sec. |5l 

We neglect all effects due to Quantum Electrodynamics (QED) on the PDF evolution. For 
reasons of consistency, QED effects (including a photon PDF, see e.g., the analysis in see [541) are 
sometimes needed in computations of cross sections including electroweak corrections at higher 
orders. Quite generally the effects are small, though. The NNLO QCD corrections to the photon's 
parton structure are known [|55ll56l and we will address this issue in a future publication. 

2.2 Power corrections 

The leading twist approximation to the QCD improved parton model is valid only at asymptotically 
large momentum transfers and the factorization underlying eqs. (12.11) and (|2.2I) is not sensitive 
to the finite hadron size effects or, equivalently, to soft hadronic scales like the nucleon mass 
A/yy. At low momentum transfer comparable to the nucleon mass such hadronic effects cannot be 
ignored and the standard factorization ansatz acquires power corrections in l/Q^. In the case of 
PDF analyses the higher twist terms are especially important for the DIS data since they cover a 
kinematical range down to ~ M^. The power corrections for the kinematics of DY data used 
in our fit (cf. Sec. [3]) are negligible due to the large momentum transfer » in this case. 
Therefore we do not consider power corrections for the DY process. 

In DIS the power corrections arise from kinematic considerations once the hadron mass effects 
are taken into account, i.e., the so-called target mass correction (TMC). The TMC can be calculated 
in a straightforward way from the leading twist PDFs within the OPE [f57l . In our analysis the 
TMC are taken into account in the form of the Georgi-Politzer prescription (ST). For relevant 
observables, i.e., the structure function F2 and the transverse one it reads 

= ^F,(^,e2)+2^ f §F2(r,e'), (2.10) 

rr Q 7 J 

and ^ 

i^2™^(^,e') = #3^2(^,e2)+6^ f ^F2(r,e2), (2.11) 

respectively, which holds up to OiM^/Q^). Here ^ = 2x/(l +y) and y = (1 +4x^Mj^/Q^)^^^ is the 
Nachtmann variable [58J. The quantities on the right hand side of eqs. (|2.10l) and (12.111) are the 
leading twist structure functions introduced in eq. (12.71) above. 

Power corrections can also arise dynamically as so-called higher twist terms from correlations 
of the partons inside the hadron. The twist-4 terms in the nucleon structure function F2 turns out to 
be non-negligible at large x [|59[I601 . Moreover, the higher twist terms in the longitudinal structure 
function appear to be necessary for the description of the NMC data at moderate x [271|6ll and the 
SLAC data on the structure function R = ailar [l62l , where ctl and ctt are the absorption cross 
sections for the longitudinally and the transversely polarized virtual photons, respectively (see also 
eq. 13.21 ). The OPE, for lepton-nucleon DIS provides the framework for the systematic classifica- 
tion of the higher twist terms referring to local composite operators of twist-4 and higher [1631 . 
Nonetheless the shapes of the higher twist terms are poorly known. Therefore they cannot be 
accounted for on the same solid theoretical footing as the leading twist contributions discussed 
in Sec. 12.11 Furthermore, both, the scaling violations and the Wilson coefficients for the various 
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higher twist contributions have not been computed to the same order in perturbation theory as for 
the leading twist part. 

Basically two strategies exist to address the issue of power corrections in the PDF analysis. The 
first one imposes kinematical cuts on the data. For DIS, these cuts are performed at high hadronic 
invariant masses = Q^{l/x- 1) + where the nucleon mass is included in the kinematical 
considerations. In this way, one aims at a data sample with reduced sensitivity to power corrections. 
Typical values for cuts on are of the order of 12 GeV^. As a drawback of this procedure one 
eliminates a rather large fraction of data at low with excellent statistical precision. A more 
serious concern, however, is due to the generally poor theoretical understanding of those non- 
perturbative QCD effects beyond leading twist factorization. One simply cannot estimate from first 
principles the region of (or W^), where power corrections can be safely neglected. Therefore, 
the present analysis (following Ii64l ) examines both, TMC and higher twist contributions, in detail 
in order to control and quantify their impact in the determination of the standard leading twist 
PDFs. 

In practice higher twist contributions are usually parameterized independently from the lead- 
ing twist one with some function of x, which is typically polynomial in x. In our analysis the 
power corrections are non-negligible for the case of the DIS data and are defined within an entirely 
phenomenologically motivated ansatz, as follows 

FYix,Q^) = F™^ix,Q^) + ^-^ + ^-^ + ..., (2.12) 

where i^™c given by eqs. (|2.10l) and (|2.1 II) . The coeflicients Hj are parameterized by a cubic 
spline with the spline nodes selected at x=0, 0.1, 0.3, 0.5, 0.7, 0.9, and 1. This choice provides 
sufficient flexibility of the coefficients Hi with respect to the data analysed and, at the same time, 
keeps a reasonable number of nodes. The values of ///(I) are fixed at zero due to kinematic 
constraints. The values of ///(O) are also put to zero in view of the fact that no clear signs of 
any power-like terms can be found in the low-x HERA data. The rest of the spline-node values 
of Hi were fitted to the data simultaneously with the PDF parameters and the value of as- We 
neglect the Q^-dependence of the higher twist operators due to the QCD evolution. Therefore 
the coeflicients Hi do not depend on Q^. This treatment could be further refined by considering 
the individual (quasi-)partonic OMEs along with their renormalization, i.e., their 2^-dependence, 
which is known for the twist-4 operators to first order in as [65J. Another complication is the 
emergence of more Bjorken-like variables Xi, the number of which grows with increasing twist. 
Experimental information on the other hand is only available for the variable x = YiiXi and Q^. We 
leave these aspects for future studies. 

With the kinematical cuts imposed on the DIS data in our analysis (cf. Sec. |3]) the twist-6 
terms are irrelevant [66] therefore the coefficients HJ"^ in eq. (12.121) are washed out. The target 
dependence of the higher twist parametrization in eq. (12.121) has been studied in [[67] - [69l . The 
isospin asymmetry in Ht is poorly constrained by the data used in our fit. It is comparable with 
zero within the uncertainties [|67l and therefore it was put to zero in our analysis. The isospin 
asymmetry in H2 is also numerically small, however, due to lower uncertainties it cannot be put 
to zero without deterioration of the fit quality. In summary, we fit three twist-4 coefficients, for 
the proton H^, for the neutron //" and for the nucleon H^, in addition to the leading twist terms. 
The impact of the power corrections on the DIS neutrino-nucleon di-muon production data used in 
our fit is marginal Il68l as well as on the inclusive charge-current data [fT8l . Therefore they are not 
considered for the case of charged-current structure functions. 
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2.3 Corrections for nuclear effects 



For our analysis we select primarily the data obtained off proton targets. However in some cases 
the necessary constraints on PDFs come only from nuclear target data. For example the neutral- 
current DIS off deuterium targets allows the separation of the large-x u- and J-quark distributions, 
which contribute to the DIS structure functions in the form of linear combination weighted with the 
quark charges. However the analysis of the deuterium data requires modeling of nuclear effects. 
They include the Fermi-motion and off-shellness of the nucleons, excess of pions in the nuclear 
matter, Glauber shadowing, etc., cf. fTOl 17111 for reviews. Among them only the Fermi-motion 
effects can be calculated with an uncertainty better or comparable to the uncertainty in the existing 
experimental data. The Fermi-motion correction is given as a convolution of the free nucleon 
structure functions with the deuterium wave function, which in turn is constrained by the low- 
energy electron-nuclei scattering data. The parameterization of the off-shell effect used in our fit 
was obtained from the analysis of the world data on DIS off heavy nuclear targets and extrapolated 
to the deuterium target. In this way we assume that the nuclear model suggested in [1701 can be 
applied to the case of light nuclei, like deuterium. This assumption has been recently confirmed 
for the case of the -^He and ^Be targets fj2\. However, in order to take a conservative estimate 
of the deuteron correction uncertainty due to the off-shellness effect we vary its magnitude by 
50%. The uncertainty obtained in this way is comparable with one given in [ItTI . Other nuclear 
effects, like shadowing and pion excesses in nuclei considered in fTOl were found numerically 
negligible for the case of deuterium. Thus our model of the nuclear effects for deuterium is based 
on the combination of the Fermi-motion and the off-shellness effects only. The nuclear corrections 
depend both on the deuterium wave function and the free nucleon structure functions, while the 
latter include the target-mass corrections and the twist-4 terms, cf. eq. (12.121 ). Due to the structure 
function dependence the value of the correction is sensitive to the fitted parameters and ideally it 
should be re-calculated iteratively in the fit. However, this approach turns out to be rather time- 
consuming. Therefore we calculate the deuteron correction once at the beginning with the PDFs 
and twist-4 terms which were obtained in [16J. Since the deuteron model employed in our fit is 
the same as one of [fT6l this approach introduces only a marginal bias into the fit. The nuclear 
correction for the representative kinematics of the deuterium data is given in Fig. |2.1| for the cases 
of the deuterium wave function obtained with the Paris potential of f73] used in our analysis and 
the Bonn potential of [74] . The difference between these two cases is marginal as compared to the 
errors in the data. The data on di-muon production in the (anti-)neutrino-nucleon scattering which 
are used in our analysis in order to constrain the strange sea distribution, were obtained on iron 
targets. In contrary to the neutral-current case a particular shape of the nuclear correction at large x 
is unimportant since data do not populate the region ofx> 0.3. At small x the neutrino-nucleon DIS 
nuclear corrections are enhanced due to the parity non-conserving part of the charged current fT5^, 
however their impact on the strange sea distribution is still smaller than its error [[68][ . Therefore 
in the modeling of the di-muon production data we employ the nuclear corrections of [[751 without 
consideration of their uncertainties. The DY data used in our fit span the moderate-x region. For 
such kinematics the nuclear effects are quite smaller than the errors in data [|76|[ and they are not 
considered in our analysis. 

3 Data 

The nucleon PDFs are usually extracted from a combination of hard-scattering data, which pro- 
vides complementary constraints on the different PDF species. A particular choice of which pro- 
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Figure 2.1: The ratio of the deuterium structure function F2 (left) and Ft (right) with account of the 
Fermi-motion and off-shellness effects of [70] calculated for the Paris potential of i73il (solid) and the Bonn 
potential of [i74i (dashes) at the momentum transfer of 20 GeV^ to the sum of those for free proton and 
neutron versus x. The shaded area around the solid line gives the uncertainty due to a variation of the 
off-shell effects by 50%. The calculations are performed at NNLO QCD accuracy using the PDFs and the 
twist-4 terms obtained in 116]. 



cesses to be used in an analysis is commonly driven by the theoretical accuracy of the data mod- 
eling and/or the experimental uncertainties in the data. In our fit we employ the data on inclusive 
DIS, the DY process, and di-muon production in neutrino-nucleon DIS. In combination they allow 
for a good separation of the quark flavors in wide range of x and provide good constraints on the 
gluon distribution at small values of x, which are mostly important for the collider phenomenology. 

3.1 Inclusive DIS 

Studies of inclusive DIS date back to the early days of QCD and since that time a wealth of the 
accurate data has been collected. The first fixed-target DIS experiments at SLAC were followed 
by data from CERN and Fermilab and then at the electron-proton collider HERA at DESY. The 
most accurate data of these experiments obtained on the proton and deuterium targets are included 
into our analysis [1311.177^841 . 

In all cases we employ the data on the inclusive cross section, which is related to the DIS 
structure functions as follows 

(f(T{x,Q^) ^ 
dxdQ^ 

with the mass of the incident charged lepton m/ and ±^^3' corresponding to different polarizations 
for the case of the charged current. The nucleon structure functions calculated 
with account of the nuclear correction described in Sec. 12.31 if relevant. Note, that higher twist 
contributions to F^^ are set to zero, cf. Sec. 12.21 In this way we provide a consistent treatment of 



xQ^ '^'-y-^- 



F\\x,Q^) + \ 



' 2mf\ 
1- 



1 
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the data, contrary to the common procedure in global PDF fits which are based on the data for the 
structure function F2. The structure functions F2 and Fj in eq. (13.11) also enter in the ratio of the 
longitudinally to transversely polarized virtual photon absorption cross sections (see e.g. [i85J), 



-1. 



(3.2) 



In order to avoid contributions from nucleon resonances and the twist-6 terms we do not include 
into the analysis any inclusive DIS data with 



< 2.5 GeV^ , 



W< 1.8GeV. 



(3.3) 



The kinematics spanned by each DIS data set used in our fit and their systematic uncertainties 
are described in the following subsections. The normalization uncertainty is a particular case of 
the systematics. However it is considered separately since very often the absolute normalization 
of the DIS experiment is not independently determined. Instead, in such cases it is usually tuned 
to a selected set of other DIS experiments, which in turn provide the absolute normalization. The 
wealth of the DIS data used in our fit allows us to extend the basis for this normalization tuning. 
Therefore, for the experiments lacking an absolute normalization we consider general normaliza- 
tion factors which are fitted simultaneously with other parameters of our data model. Within this 
approach we introduce free normalization parameters for the separate early SLAC experiments 
of [|77] - l80l and for the NMC data of [84J at each beam energy. The errors in the normalization 
factors obtained in our fit are included into the general covariance matrix calculation. In this way 
we account for the impact of the absolute normalization uncertainty in the data on the PDFs, the 
higher twist terms and on the value of a^. 

Our procedure for the treatment of the DIS data normalization in PDF fits differs substantially 
from other approaches. For instance, in the MSTW PDF fit [22 J free normalization parameters 
are introduced for all data sets, including even those where the absolute normalization has been 
determined experimentally. Other PDF fits also commonly employ the NMC data averaged over 
the beam energies and combined data from the SLAC experiments, rather than the respective indi- 
vidual data sets. 

3.1.1 HERA 

In our analysis we use the HERA data on the inclusive neutral-current and the charged-current cross 
sections [18|. This sample was obtained by a combination of the run I data of the HI and ZEUS 
experiments, and includes in particular the data of [f86ll87l used earlier in the ABKM09 fit [fT6l . 
The HERA data span the region of up to 30000 GeV^. However, we impose an additional 
cut of < 1000 GeV^ on the neutral-current sample. This allows to avoid the complexity due 
to the Z-boson exchange contribution. At the same time, the high-Q^ part of those data displays 
only a poor sensitivity to our PDF fit. The normalization uncertainty in the HERA data of [18] is 
0.5%, much better than one in the HERA data of [86, 87J. In particular due to the improvement 
in normalization the new data of [18J somewhat overshoot the previous HI data of [86J. The 
total number of correlated systematic uncertainties in the HERA data of [[T8ll is 1 14, including the 
uncertainties due to the combination procedure and the general normalization. Many of them are 
improved as compared to the separate experiment samples as a results of cross-calibration in the 
process of combination. 
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A complementary set of the inclusive HERA data was obtained by the HI collaboration in the 
run with a reduced collision energy [31 J. These data are particularly sensitive to the structure func- 
tion Fl and thereby to the small-x shape of the gluon distribution. The normalization uncertainty 
in the low-energy data of [[3T|| is 3%. The point-to-point correlated systematic uncertainties come 
from 8 independent sources and there is also a number of uncorrelated systematic uncertainties in 
the data. 

3.1.2 BCDMS 

The BCDMS data of [|82l[83l used in our fit were collected at the CERN muon beam at energies 
of 100, 120, 200, and 280 GeV for the incident muons. Due to the use of both proton and deu- 
terium targets in the same experiment these data facilitate flavor separation of PDFs at large x. The 
BCDMS absolute normalization was monitored for the beam energy of 200 GeV. The general nor- 
malization uncertainty in the data due to this monitor is as big as 3%. The absolute normalization of 
the data obtained at the beam energies of 100, 120, and 280 GeV was calibrated with respect to the 
case of the beam energy setting of 200 GeV. The additional normalization uncertainty due to this 
calibration ranges from 1% to 1.5% depending on the beam energy. Other systematic uncertainties 
in the BCDMS data stem from 5 uncorrelated sources with the most important contributions due 
to incident and scattered muon energy calibration and the spectrometer resolution. 

3.1.3 NMC 

The NMC experiment was performed like BCDMS at the CERN muon beam at incident muon 
energies of 90, 120, 200, and 280 GeV. However, the NMC data span lower values of x and as 
compared to BCDMS and overlap with the HERA data at the edge of the respective kinematics. 
We use in our fit the NMC cross section data of [84] for the proton and deuterium targets. Due 
to better coverage of the small-.x; region those data are also sensitive to the isospin asymmetry in 
the sea distribution. The absolute normalization for the NMC data of [84] was determined from 
tuning for each particular energy setting separately to the BCDMS and SLAC data, which overlap 
partially with NMC. This tuning in [|84|| was based on an empirical data model motivated basically 
by leading-order QCD calculations. 

In our analysis, therefore, we fit the NMC normalization factors for each incident beam energy 
and target simultaneously with the other parameters. In this manner, we ensure consistency with 
our data model, which in particular includes QCD corrections up to the NNLO, see [f27| for a 
detailed study of the impact of the NNLO QCD corrections on the interpretation of the NMC data. 
The normalization factors obtained in the NNLO variant of our fit are given in Tab. 13.11 In general 
they are within the uncertainty of 2% quoted for the NMC data in [ 88 1 . However, the normalization 
factors for the proton target are somewhat larger than for the case of deuterium. This is explained 
by impact of the HERA data of [[T8ll . which slightly overshoot the NMC data in the region of 
their overlap. The systematic uncertainties in the NMC data are due to the incident and scattered 
muon energy calibration, the reconstruction efficiency, acceptance, and the electroweak radiative 
corrections. Some of the systematic uncertainties are correlated for all data, some of them between 
the proton and deuterium data, and some between beam energies (cf. [ l84l | for details). In summary 
this gives 12 independent sources of systematic uncertainties for the NMC data used in our fit. 
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Beam energy (GeV) 


proton 


deuterium 


90 


1.012(12) 


0.990(12) 


120 


1.026(11) 


1.005(11) 


200 


1.034(12) 


1.014(11) 


280 


1.026(11) 


1.007(11) 



Table 3.1: The NMC normalization factors obtained in our NNLO fit for different incident beam energies 
and targets. 



3.1.4 SLAC 

The SLAC experiments used in our fit and the number of data points for each experiment after 
the cut of eq. (|3.3I) are listed in Tab. 13. 2[ The last and most elaborated in this series is experiment 
E-140 n8ll . In particular, it took advantages of the improved electroweak radiative corrections 
and the accurate determination of the data absolute normalization, which is as big as 1.8% for the 
deuterium sample. Other point-to-point correlated systematic uncertainties are due to background 
contamination, the spectrometer acceptance, and the electroweak radiative corrections [|89ll . The 
rest of systematic error sources for the experiment E-140 are uncorrelated. 



Experiment 


Target 


NDP 


NSE 


Normalization 
(our fit) 


Normalization 
(Ref. il) 


E-49a|77| 


proton 
deuterium 


59 
59 


3 
3 


1.022(11) 
0.999(10) 


1.012 
1.001 


E-49b [771 


proton 
deuterium 


154 
145 


3 
3 


1.028(10) 
1.008(10) 


0.981 
0.981 


E-87 [77] 


proton 
deuterium 


109 
109 


3 
3 


1.032(10) 
1.017(10) 


0.982 
0.986 


E-89a DTI 


proton 
deuterium 


77 
71 


4 
5 


1. 
1. 


0.989 
0.985 


E-89b |79| 


proton 
deuterium 


90 
72 


3 
3 


1.016(10) 
0.996(10) 


0.953 
0.949 


E-139 [80J 


deuterium 


17 


3 


1.014(10) 


1.008 


E-140 [f8T1l 


deuterium 


26 


5 


1. 


1. 



Table 3.2: The list of SLAC experiments used in our fit (first column: the experiment number; third 
column: the number of data points (NDP) used in the fit; fourth column: the number of correlated systematic 
errors (NSE) in the data; fifth column: the normalization factor applied to the data in our fit; sixth column: 
the normalization factor applied to the data in the re-analysis of I62J ). Note, that the normalization factors 
of the fifth column apply to the data, which were re-normalized in ll62i by the factors given in the sixth 
column. 

The earlier SLAC data used in our fit are collected with various experimental setups and data 
processing chains. In particular, the electroweak radiative corrections applied to the data differ 
in details, various methods are used to determine absolute normalization of the data, etc. To 
overcome this diversity the early SLAC data of iTTl - ISOll were reanalysed within a uniform approach 
and the leveled set of the SLAC data was obtained in [|62l| . As a part of this leveling the absolute 
normalization factors for the data of [l77] - [80ll were calibrated with the help of the E- 140 data of ifSTI . 
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Due to the lack of the E-140 proton data this calibration is straightforward only for the deuterium 
case. The proton data normalization tuning was performed in two steps. First, the normalization of 
the deuterium data of experiment E-49b was determined with the help of the E-140 deuterium data. 
Then, the proton data normalization for all other experiments was tuned to the E-49b proton data, 
assuming equal normalization for the proton and deuterium samples of the E-49b experiment. The 
data of experiment E-89a are kinematically separated from other SLAC experiments considered. 
Therefore, their normalization tuning was based on the elastic scattering samples obtained in the 
experiments E-89a, E-89b, and E-140 (cf. [62] for the details). In view of the fact, that we do 
not include elastic data in the fit we keep the normalization of the SLAC experiment E-89a at the 
value obtained in [|62ll . At the same time in order to take into account the uncertainties in the E- 
89a data normalization we add to those data the general normalization uncertainty of 2.8% and an 
additional normalization uncertainty of 0.5% for the case of deuterium, which are quoted in ||62l . 
The normalization factors for the early SLAC experiments of fTTl 179118011 are considered as free 
parameters of the fit. The SLAC normalization factors for the NNLO variant of our fit are given 
in Tab. 13.21 in comparison with the ones of [[62ll , which were obtained with an empirical QCD- 
motivated model of the data. The deuterium normalization factors obtained in our fit are in a good 
agreement with the ones of Il62|| . For the proton target case our normalization factors are somewhat 
bigger, in particular due to wider set of data is used for the normalization tuning in our case. 

3.2 Drell-Yan process 

The data for the DY process provide a complementary constraint on the PDFs. In particular they 
allow to separate the sea and the valence quark distributions in combination with the DIS data. 
We use for this purpose the data obtained by the fixed-target Fermilab experiments E-605 [.90J and 
E-866 jm. 

The experiment E-605 collected proton-copper collisions data at the center-of-mass energy of 
38.8 GeV for di-muon invariant masses in the range of 7 17 GeV. At this kinematics the DY 
data are sensitive to the PDFs down to x ~ 0.03. The normalization uncertainty in the E-605 data 
is 15%. However other systematic uncertainties in the data are not fully documented in f90'|. The 
point-to-point correlated systematic is estimated as -i-10% for low di-muon masses and -10% for 
higher masses. Due to lacking details in ||90l we assume a linear dependence of this systematic 
error on the di-muon mass. Additional uncorrelated systematic uncertainties in the E-605 data due 
to the Monte Carlo acceptance calculation are combined with the statistical ones in quadrature. 

The data of the E-866 experiment on the ratio of the proton-proton and proton-deuterium col- 
lision cross sections ||9T| are particularly sensitive to the isospin asymmetry of the sea quark dis- 
tributions. The absolute normalization uncertainty cancels in this ratio. Other E-866 systematic 
uncertainties stem from 5 independent sources with the biggest contributions due to the deuterium 
composition and the event detection and reconstruction. The unpublished data on the absolute DY 
cross sections for the proton and deuterium targets are also available [92] . However these data are 
in poor agreement with the DIS data (cf. [93] for a detailed comparison). Therefore, we do not 
employ the data of (921 in our fit. Note, that in the MSTW PDF fit JH the E-866 data on the 
absolute cross sections are shifted upwards by 8.7% in order to bring them into agreement with the 
other data sets. 
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3.3 Di-muon production in vN DIS 

The production of di-muons in neutrino-nucleon collisions provides unique information about the 
strange sea distribution in the nucleon. One of the muons produced in this reaction may be resulting 
from the decay of a charmed hadron. Thus, the production of the c-quarks in neutrino-nucleon 
collisions is directly related to initial-state strange quarks. Therfore, by relying upon a c-quark 
fragmentation model one can determine the (anti-)strange-sea distribution from the data on di- 
muon production in an (anti-)neutrino beam. The details of the fragmentation model are quite 
important in this context due to kinematic cuts imposed to suppress a background of muons coming 
from the light mesons. Herewith the absolute normalization of the model is defined by the semi- 
leptonic branching ratio 5^ of the charmed hadrons. The value of 5^ is poorly known due to the 
uncertainty in the hadronic charm production rate for the neutrino-nucleon interactions. On the 
other hand, the value of 5^ is also constrained by the di-muon data themselves [68,94]. Therefore, 
for consistency, we fit the value of 5^ simultaneously with the PDF parameters imposing available 
independent constraints on 5^ coming from emulsion experiments (cf. [[68^1 for the details). 

We use in the fit the di-muon data provided by two Fermilab experiments, CCFR and NuTeV, 
and corrected for the cut of 5 GeV imposed on the muon decay energy in order to suppress the 
light-meson background [|95ll96l . The data of the NuTeV experiment were normalized through 
the use of the inclusive single muon event rates. Therefore, the normalization error in the data 
is marginal and it is not considered in our fit. Besides, 8 independent sources contribute to the 
point-to-point correlated systematic uncertainties. The neighboring NuTeV data points are also 
correlated due to smearing of the kinematic variables. These correlations are not documented 
in L95J. Instead the errors in the data are inflated in such a way that the fit of a model to the data 
with inflated errors is equivalent to the regular fit with account of the data correlations (cf. Il95l for 
details). The average data error inflation factor is about 1.4. Therefore, the normal value of for 
the inflated-error fit is about one half of the NDP. The CCFR data of f95l were processed similarly 
to the NuTeV ones. In particular, the errors in the data were also inflated by factor of about 1.4 
in order to take into account the data point correlations. However, only the combined systematic 
errors in the CCFR data are available. In view of the lack of any detailed information about the 
systematic error correlations we employ in our fit the combined systematic errors quoted in [|95l 
assuming them to be fully point-to-point correlated. 

4 Results 

We are now in a position to present the results of our analysis ABMl 1 to NLO and NNLO in QCD 
fortif = 3 in a FFNS. The PDF sets for rif = 4 and 5 are then generated by matching as described 
above and we will comment on the changes in the PDFs obtained compared to the ABKM09 
set [fT6l . In addition to the fit results and the co variance matrix for the correlations of the fit 
parameters, we also present the pulls for separate experiments, which reflect the compatibility of 
these data sets with respect to the experimental systematics. The discussion of the value of the 
strong coupling obtained in ABMll is supplemented by a compilation of as{Mz) determinations 
in NLO and NNLO analyses extending [97 1. For a valence distribution we also compute the lowest 
Mellin moment of our PDFs and compare with the latest available data from lattice simulations. 
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4.1 PDF parameters 



In the new analysis the shape of the PDFs has been updated and the number of fit parameters 
has been slightly enlarged compared to ABKM09 lfT6l . In detail, we are using the following 
parametrizations at the starting scale yu^ = 2^ = 9.0 GeV^ in the scheme with «/ = 3 flavors, 



•^'?v(^, Qq) 

xus(x, qI) = xusix, qI) 
xA(x, qI) = xdsix, qI) - xusix, qI) 
xsix, qI) = xsix, qI) 



Aajc"*(1-jc)*Ax^^^^\ 



= A,x"-'i\-xf 



xgix,Q'o) = Agx"^^(l-xr^x' 



(4.1) 

(4.2) 
(4.3) 
(4.4) 
(4.5) 



where q = u,d and 6qq' denotes the Kronecker function in eq. (|4.1I) and the strange-quark distribu- 
tion is taken to be charge-symmetric, cf. [168]| . The polynomials P(x) in eqs. (I4.1I) - (I4.5I) are given 
by 



Pqv(x) 



Pus(x) 
PAix) 

PJx) 



7hqX + 72,qX^ +y3,qX^ , 

\nx{\ +n,usX + y2,usX^), 

7l,AX, 
71,8^- 



(4.6) 
(4.7) 
(4.8) 
(4.9) 



The new functional form with the additional parameters 73 „ and y^^us provides sufficient flexibility 
in the small-.x; M-quark distribution with respect to the analysed data and we have checked that no 
additional terms are required to improve the quality of the fit. All 24 PDF parameters are given in 
Tab. 14. II together with their Icr uncertainties computed from the propagation of the statistical and 
systematic errors in the data, cf. App. |Al Note that the normalization parameters for the valence 
quarks, A^^, and gluons, Ag, are related to the other PDF parameters due to conservation of fermion 
number and of momentum, respectively. 





a 


b 


ri 


72 


73 


A 


My 


0.712 ±0.081 


3.637 ±0.138 


0.593 ± 0.774 


-3.607 ± 0.762 


3.718 ± 1.148 




dv 


0.741 ±0.157 


5.123 ±0.394 


1.122 ± 1.232 


-2.984 ± 1.077 






Us 


-0.363 ± 0.035 


7.861 ± 0.433 


4.339 ± 1.790 




0.0280 ± 0.0036 


0.0808 ± 0.0122 


A 


0.70 ± 0.28 


11.75 ± 1.97 


-2.57 ±3.12 






0.316 ±0.385 


s 


-0.240 ± 0.055 


7.98 ± 0.65 








0.085 ± 0.017 


8 


-0.170 ±0.012 


10.71 ± 1.43 


4.00 ±4.21 









Table 4. 1 : The parameters of the PDFs in eqs. (I4.1l) - (l4.5b and their Icr errors obtained in the scheme with 
n/ = 3 flavors. 



As in our previous analysis ABKM09 [fT6ll . the small-.x; exponent for the difference between 
the up- and the down-quark sea is fixed to a/s. = 0.7 in eq. (14.31) as an ansatz, because of lacking 
neutron-target data in this region of small values of x. This is in agreement with the values ob- 
tained for the small-.x: exponents of the valence quark distributions. The uncertainty on is then 
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determined with help of an additional pseudo-measurement of aA = 0.7 + 0.3 added to the data set 
(and with the error on released) in order to quantify the impact on the other parameters of the 
PDF fit. This provides us with the result given in Tab. 14.11 The value of the charmed-hadron semi- 
leptonic branching ratio obtained from our NNLO fit is 0.0917 + 0.0034. This is in agreement 
with the earlier determination ||68l within the errors. 
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Figure 4.1: The pulls versus momentum transfer for the HERA neutral-current inclusive DIS cross 
section data of 1 18 1 binned in x with respect to our NNLO fit. The data points with different inelasticity y 
still may overlap in the plot. The inner bars show statistical errors in data and the outer bars the statistical 
and systematic errors combined in quadrature. 

The three other parameters of our fit to be discussed in detail in Secs. l4.2l and l4.3[ are the strong 
coupling constant Us in the MS -scheme and the heavy-quark masses and mi,, which we take in 
the MS -scheme as well. The latter represents a novel feature of our analysis as all previous PDF 
determinations have always used the pole mass definition for nic and m^. As an advantage, we can 
constrain the central values of both, mdmc) and mb(mb) directly to their particle data group (PDG) 
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Figure 4.2: The same as Fig. 14. II for the pulls of the HERA charged-current inclusive DIS cross section 
data of [181 binned in the momentum transfer in units of GeV^ versus x (squares: positron beam; circles: 
electron beam). 
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Figure 4.3: The same as Fig. 14. II for the pulls of the HI neutral-current inclusive DIS cross section data 
of ||3T| binned in the momentum transfer in units of GeV^ versus x. 
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results [W| without having to rely on a perturbative scheme transformation between a running MS - 
and a pole mass. It is well known that at low scales such as fi ^ nic this scheme transformation 
is poorly convergent in perturbation theory. Thus, in the present analysis, we add the following 
pseudo-data as input 

mdmc) = 1.27 + 0.08 GeV, = 4.19 + 0.13 GeV, (4.10) 

and, subsequently, release the uncertainty of the quark masses to test its sensitivity to the other 
PDF parameters. The value for as{Mz) on the other hand is determined entirely from data in the 
fit, cf. Sec. 14. 2[ The 24 PDF parameters of Tab. 14. 1[ Us, nic and provide us in our analysis in 
total with 27 correlated parameters. Their covariance matrix is presented in Tabs. lA[l"l - IA[3[ 
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Figure 4.4: The same as Fig. [JT] for the pulls of the BCDMS inclusive DIS cross section data of JH 
for the proton target (left) and for the deuterium target (right). 



It is instructive to study the pulls of the individual data sets included in the fit. This provides a 
mean of assessing the quality of the fit in detail and allows for an investigation of specific kinemat- 
ical regions. In Figs. 14.11 - 14. 31 we display the detailed dependence of the pulls on the momentum 
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Figure 4.5: The same as Fig. |4?T]for the pulls of the NMC inclusive DIS cross section data of lIMl for the 
proton target (left) and for the deuterium target (right). 



transfer and x for the HERA NC and CC inclusive DIS cross section data of IfTSl as well as the 
low data of [31 1 with respect to our NNLO fit. We find overall a very good description of the 
data, even at the edges of the kinematical region of HERA, i.e., at smallest values of x and largest 
values of Q^. The respective values for the fit at NLO and NNLO are given in Tab. 14.21 

Next, in Figs, g^ffl we show the respective pulls of the BCDMS [l82l[83l, NMC MM and 
SLAC |[77] - [8Tll inclusive DIS cross section data as a function of x and binned in the momentum 
transfer Q^. Again, our fit provides a very good description (see Tab. 14.21 for the values), 
especially at low thanks the phenomenological ansatz for the structure functions with the higher 
twist terms of Tab. I4.3[ 

In Fig. 14.71 we plot the data for the (anti-)neutrino induced di-muon production cross section of 
||95l which constrains the strange PDF. We give both, the pulls for the NuTeV and for the CCFR 
experiment. Finally, in Fig. 14.81 we display the DY cross section data of [90l|91]| which depends 
on the muon pair rapidity Y and the invariant mass M^^ of the muon pair and which assists in the 
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Figure 4.6: The same as Fig. l4.1l for the pulls of the SLAC inclusive DIS cross section data of l77l - [8T]| for 
the proton target (left) and for the deuterium target (right). 



flavor separation of the PDF fit. It is obvious from Figs. |4.71l4.8l and Tab. |4.2| that we achieve again 
a very good description in all cases. 

The last missing piece of information on the PDF fit concerns the shape of the higher twist 
terms for the inclusive DIS structure functions introduced in eq. (|2.12l) . As outlined in Sec. 12.21 
we fit three twist-4 coefficients for a complete description of both, proton and nucleon targets. In 
detail, the proton H^, the non- singlet H^^ = H^- and the proton Hj contribute 15 parameters 
in total and we assume //"^ = 0. The respective coefficients are listed in Tab. 14.31 and shown 
in Fig. |4.9[ all in units of GeV"^, They are in agreement with the earlier results of [^9l[60ll up 
to the parameterization of the HT contribution (compare eq. (12.121) with eq. (35) of [601). The 
magnitude of the HT terms reduces from the NLO to the NNLO case, see Fig. 14. 9[ however the 
change is comparable with the coefficient uncertainties and the NNLO twist-4 coefficients do not 
still vanish, in line with the results of ||60l . The non-singlet twist-4 term in Ft is comparable to 
zero within uncertainties therefore it was fixed at zero in our analysis as discussed in Sec. 12.21 The 
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Figure 4.7: The same as Fig. 14. II for the pulls of neutrino (left) and anti-neutrino (right) induced di-muon 
production cross section data of Ii95il (circles: NuTeV experiment, squares: CCFR experiment). 
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non-singlet twist-4 term in F2 is negative at x < 0.5. This is also in line with earlier results ll59l . 
again taking into account the difference in the HT parameterizations. The HT terms are mostly 
important at small hadronic invariant mass W. This was confirmed in a comparison of the low-V7 
JLAB data HIOOH with predictions based on the ABKM09 PDFs with account of the twist-4 terms, 
which were extracted in the analysis of [16J similarly to the present one. However even with the 
cut of > 12.5 GeV^ as commonly imposed in global PDF fits the HT terms are numerically 
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SLAC-E-139 iOl 


17 


11 


11 




SLAC-E-140 [81] 
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28 


29 


Drell-Yan 


FNAL-E-605 [,90il 
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FNAL-E-866 JUl 
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52 


55 


DIS di-muon 


NuTeV [95! 


89 


46 
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CCFR [95J 


89 


61 


62 


Total 




3036 


3391 


3378 



Table 4.2: The value ofx^ obtained in the NNLO and NLO fits for different data sets. 




Figure 4.9: The central values (solid line) and the Icr bands (shaded area) for the coefficients of the twist-4 
terms of eq. ( 12.121 ) in the inclusive DIS structure functions obtained from our NNLO fit (left panel: Fi of 
the proton, central panel: Fj of the proton, right panel: non-singlet f 2)- The central values of the twist-4 
coefficients obtained from our NLO fit are shown for comparison (dashes). 



important for the region of x < 0.3, which is not affected by this cut. 

From Fig. 14.101 it is evident, that calculations, which are based on our NNLO PDFs, but do 
not include the HT terms, systematically overshoot the SLAC data at x < 0.3 due to the HT terms 
being negative in this region. The value of;^^^/NDP for the SLAC data at 

> 12.5 GeV^ > 2.5 GeV^, (4.11) 

is 699/246 in this case. This is much worse than the value of;t'^/NDP = 292/246 obtained in our 
fit for the same subset of data. We have also performed similar comparisons taking the published 
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= 0.7 
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-0.029 + 0.006 
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= 0.9 


0.0071 +0.0026 


0.0009 + 0.0041 


0.0002 + 0.0015 



Table 4.3: The paiameters of the twist-4 contribution to the DIS structure functions in eq. (12.121 ) for the fit 
to NNLO accuracy in QCD. 



3-flavor NNLO MSTW ^ and NLO NN21 B PDFs as an input of our fitting code. The NNLO 
MSTW and NLO NN21 predictions obtained in this way without accounting for the HT terms and 
with the cut of eq. (14.1 II) imposed lead to poor a description of the SLAC data, cf. Fig. 14.101 For 
the case of NN21 PDFs the agreement with the data is particularly bad, with an off-set reaching 
up to ~ 10% atx~ 0.15 and a value ofx^/^^P = 583/246. The MSTW predictions are in better 
agreement with the data, however the value of x^/NDF = 514/246 is also far from ideal in this 
case, obviously due to the missing HT terms. Note, that for the comparison performed without the 
twist-4 terms the ABM 11 value ofx^ is worse than ones of MSTW and NN21 since those PDFs 
are obtained disregarding the HT terms. This also shows that parts of the twist-4 terms obtained in 
our fit are effectively absorbed into the MSTW and NN21 PDFs. 

The relative contribution of the higher twist terms to the ratio R = (TlIctt in eq. (13.21) is partic- 
ularly important reaching up to one half at moderate x HIOIL The value of R calculated including 
the NNLO QCD corrections and the twist-4 terms of Tab. I4.3l is in reasonable agreement with the 
SLAC data on R [62J and the parameterization of those data R\99q, see cf. Fig. 14.1 II The latter is 
based on the empirical combination of the QCD-like terms with the twist-4 and twist-6 terms, pre- 
tending to describe the data down to scales 2^ ~ 1 GeV^. Due to the twist-6 term, which provides 
saturation of i?i99o at small Q, the shape of i?i99o is somewhat different from our calculation, while 
both agree with the data at 2^ > 2 GeV^ within the errors. The leading twist NNLO contribution 
to R undershoots the full calculation by a factor of 1.5 - 2, depending on x, cf. Fig. 14.1 II The 
leading-twist NLO calculations based on the 3-flavor NN21 PDFs are in a good agreement with 
our NNLO leading-twist term and go by factor of 1 .5 - 2 lower than the data as well. The leading- 
twist NNLO calculations for the 3-flavor MSTW PDFs are larger than the NLO NN21 ones and 
are in better agreement with the data. Note, that this is related to the fact that the data on R of [62]| 
are included into the MSTW fit allowing for a better description of the SLAC cross section data as 
compared to the NN21 case, see Fig. 14. 101 and the related discussion above. At the same time this 
leads to an effective absorption of the twist-4 terms into the fitted PDFs. 

The NNLO ABMl 1 PDFs obtained in the present analysis are compared in Fig. 14.121 with our 
earlier ABKM09 PDFs. The biggest change between these two sets is observed for the small-x 
gluon and sea distributions. Firstly, this change happens since the HERA NC inclusive data [13| 
used in the present analysis lie by several percent higher than the HERA data of [86.87] used in the 
ABKM09 analysis, due to improvements in the monitor calibration. Secondly, the small-x PDFs 
are particularly sensitive to the treatment of the heavy-quark electro-production and, therefore, they 
change due to the NNLO corrections and the running-mass scheme implemented in the present 
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Figure 4.10: The same as Fig. 14.61 without the HT terms taken into account and for various 3-flavor PDFs 
(left panel: present analysis, right panels: MSTW (TT]). The NLO calculations based on the 3-flavor NLO 
NN21 PDFs [99] are given for comparison (central panels). Only the data surviving after the cut of eq. (14.1 II ) 
are shown; the proton and deuterium data points are superimposed. 



PDF fit. Other ABMl 1 PDFs are in agreement with the ABKM09 ones within the uncertainties. 

The NNLO PDFs obtained by other groups are compared with the NNLO ABMll PDFs in 
Fig. 14.131 The agreement between the various PDFs is not ideal, a fact that may be explained by 
the differences in the data sets used to constrain the PDFs, by the factorization scheme employed, 
by the treatment of the data error correlation and so on. A detailed clarification of these issues is 
beyond the scope of the present paper. Therefore we discuss only the most significant differences, 
e.g., the gluon distributions at small x < 0.001, which are quite different for all PDF sets considered 
in Fig.l4l3l 

To that end, we compare in Fig. 14. 141 the small-x data on obtained by the HI collabora- 
tion IITlJ with the predictions based on these PDFs. The Fi data are quite sensitive to the small-x 
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Figure 4. 1 1 : The shaded area gives Icr band of the NNLO predictions for the ratio R = ctl/o-t, cf. eq. (I3.2I ). 
based on the ABM 1 1 PDFs and the twist-4 terms obtained from our fit at different values of x versus the 
momentum transfer Q^. The central values of the NNLO predictions for R based on the MSTW PDFs Il22l 
(dashed dots), the NNLO predictions based on the ABM 11 PDFs (solid line), and the NLO predictions 
based on the NN21 PDFs [99] (dots), all taken in the 3-flavor scheme and without twist-4 terms, are given 
for comparison. The data points show values of R extracted from the SLAC proton and deuterium data l62l 
with the empirical parameterization of those data /?i99o obtained in [621 superimposed (dashes). 

gluon PDFs. Moreover, in order to provide a consistent comparison all predictions are taken in 
the running-mass 3-flavor scheme with the heavy-quark masses of eq. (14.101) and with the 3-flavor 
PDFs. The HI data on are in a good agreement with the NNLO ABMl 1 predictions. Although 
these data were not included into the fit of [20] they are also in a good agreement with the NNLO 
JR09 predictions. The NNLO MSTW and NLO NN21 predictions on the other hand miss the HI 
data. Thus, the latter can be used to consolidate the small-;c behavior of the gluon PDFs provided 
by different groups. Likewise, the SLAC DIS cross section data of llTTl - fSTl can also be of help 
in consolidating the results of the different PDF fits. As one can see in Fig. 14.101 the MSTW and 
NN21 predictions systematically overshoot the SLAC data at x ~ 0.2. As we discussed above, this 
happens due to the omission of the higher twist terms. Once the latter are neglected in the fit, the 
power corrections are partially absorbed in the leading twist PDFs. Therefore, this discrepancy 
is evidently also related to the difference of those PDFs with the ABMl 1 ones at moderate x, cf. 
Fig. 14. 131 On the other hand, the ABMl 1 large-x gluon distribution goes lower than the NN21 and 
MSTW ones, because we do not include the Tevatron inclusive jet data into the fit, cf. IIT02i and 
Sec. [541 

Another striking difference in Fig. l4.13l is related to the strange sea distribution, which is com- 
monly constrained by the data on the di-muon production in the vA'^ DIS in all PDF fits considered. 
Nonetheless, for the NN21 and JR09 sets it goes significantly lower at .x; > 0.02 than for the MSTW 
and ABMl 1 ones. The difference between the NN21 and ABMl 1 strange sea distributions appears 
to be due to eq. (34) of [99| for the di-muon production cross section, which contains an additional 
factor of (1 +ml/Q^) as compared, e.g., to eq. (3) of [|38l employed in our analysis, cf. also [|50l . 
At small this factor reaches a numerical value of 2 and the strange sea is suppressed corre- 
spondingly in the fit to the data. On the other hand, the discrepancy between the results of JR09 
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Figure 4.12: The Icr band for the 4-flavor NNLO ABKM09 PDFs im at the scale of = 2 GeV versus x 
(shaded area) compared with the central values for ones of this analysis (solid lines). 



and ABMl 1 in Fig. 14. 131 can directly be traced back to the ansatz s{x, Q^) = s{x, Qq) = of JR09, 
i.e., the assumption of vanishing strangeness at the low starting scale of 2^ < 1 GeV^ in the dy- 
namical valence-like PDF model of JR09. This is different from ours, cf. eq. (14.41) . Moreover, 
JR09 has not used the data on di-muon production in neutrino-nucleon collisions in their fit, see 
Sec. [331 

Finally, the difference in the large-x non-strange quark distributions appears partly due to the 
general normalization of the data, which is often a matter of choice in the PDF fits. In order to 
quantify impact of the choice of the data normalization on the PDFs (cf. Sec. [3]), we have performed 
a variant of our NNLO fit with the same normalization factor settings as employed in the MSTW 
fit [|22|. The relative difference between the u- and J-quark distributions obtained in this variant of 
the fit and our nominal one is displayed in Fig. 14. 151 Clearly, the impact of the data normalization 
choice is most pronounced at large x, where the trend is different for the cases of d- and M-quarks. 
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Figure 4.13: The 1 a band for the 4-flavor NNLO ABM 1 1 PDFs at the scale of - 2 GeV versus x obtained 
in this analysis (shaded area) compared with the ones obtained by other groups (soUd lines: JR09 |[20ll . 
dashed dots: MSTW [22], dashes: NN21 |i99J). 

Therefore the effect is amplified in the ratio d/u which is important for the interpretation of the 
charged-lepton asymmetry data from hadron colliders, cf. Sec. 15.11 As shown in Fig. I4.15[ the 
relative difference in the ratio d/u reaches up to 5% at .x: ~ 0.5. 

4.2 Strong coupling constant 

For a precision determination of as(Mz) the fit of systematically compatible data sets is a necessary 
prerequisite. Here the individual data sets determine the average in such a way that their individual 
effect is closely compatible with the central value within the errors. Enhancing the precision from 
NLO to NNLO, and in some cases to even higher orders, the central values and the values obtained 
for the individual data sets both stabilize. Moreover, the values of as(Mz) obtained by individual 
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Figure 4.15: The Icr band for the 3-flavor NNLO w-quark distribution (left panel), J-quark distribution 
(central panel), and the diu ratio (right panel) at the scale ofji-l GeV versus x in comparison to the central 
value of the fit variant with the data normalization changed from the settings of Sec. [3] to the ones of Ii22l 
(solid curves). 



experiments, capable to measure as from their data alone, have to be consistently reproduced. 
One observes a decreasing sequence of differences |Aq;^(Mz)| between the sequential orders, cf. 

e.g., miM. 
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Figure 4. 16: The -profile eq. (|AI1I) as a function of OsiMz) in the present analysis. At NLO (circles) and 
NNLO (squares). 



Experiment 


a,(Mz) 


NLO,.,p 


NLO 


NNLO 


BCDMS 

NMC 

SLAC 

HERA comb. 
DY 


0.1111+0.0018 

117 +0.011 

-0.016 


0.1150 + 0.0012 
0.1182 + 0.0007 
0.1173+0.0003 
0.1174 + 0.0003 
0.108 +0.010 


0.1084 + 0.0013 
0.1152 + 0.0007 
0.1128 + 0.0003 
0.1126 + 0.0002 
0.101 +0.025 


ABMll 




0.1180 + 0.0012 


0.1134 + 0.0011 



Table 4.4: Comparison of the values of a,{Mz) obtained by BCDMS US and NMC |[T03l at NLO with the 
individual results of the fit in the present analysis at NLO and NNLO for the HERA data HSHSD, the NMC 
data m, the BCDMS data ElllMl, the SLAC data Il62ll77ti8n. and the DY data l[90ll9T]l . 



Experiment 






NLO 


NNLO* 


DO 1 jet 
DO 2 jet 

CDF 1 jet (cone) 
CDF 1 jet (k^) 


1161+™! 


0.1190 + 0.0011 
0.1174 + 0.0009 
0.1181+0.0009 
0.1181+0.0010 


0.1149 + 0.0012 
0.1145 + 0.0009 
0.1134 + 0.0009 
0.1143 + 0.0009 


ABMll 




0.1180 + 0.0012 


0.1134 + 0.0011 



Table 4.5: Comparison of the values of a^iMz) obtained by DO in II 1041 with the ones based on including 
individual data sets of Tevatron jet data II105I4108II into the analysis at NLO. The NNLO* fit refers to the 
NNLO analysis of the DIS and DY data together with the NLO and soft gluon resummation corrections 
(next-to-leading logarithmic accuracy) for the 1 jet inclusive data, cf. Ill02i|109l . 

In the present analysis based on the measured scattering cross sections and with account of 



■ NNLO 
• NLO 
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Figure 4.17: The -profile versus the value of a^iMz), eq. (lAlll) . for the data sets of Tab. l44l all obtained 
in variants of the present analysis with the value of fixed and all other parameters fitted (soUd lines: 
NNLO fit, dashes: NLO fit). 

higher twist contributions, c.f. also [,27il we obtain from the data sets described in Sec. [31 

as(Mz) = 0.1180 + 0.0012 atNLO, (4.12) 

a,(Mz) = 0.1134 + 0.0011 at NNLO. (4.13) 

The value at NNLO is shifted by AasiMz) = 0.0046 downward if compared to the NLO value. This 
range of uncertainty is well compatible with the scale uncertainty observed in a variation of the 
factorization and renormalization scales at NLO of 0(0.0050), cf. [|111| . The present data allow 
a measurement of OsiMz) with an accuracy of 0(1%). Therefore NNLO analyses are manda- 
tory, since the NLO results exhibit much too large theory errors. The response to the fitted a^- 
dependence is measured using the;t'^ functional of eq. (lAlll) . In Fig. l4.16l the dependence ofx^ of 
as{Mz) is illustrated at NLO and NNLO. 
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Figure 4. 18: The -profile versus the value of OsiMz), eq. (|A[lb . for the data sets of Tab. l44l all calculated 
with the PDF and HT parameters fixed at the values obtained from the fits with asiMz) released (sohd lines: 
NNLO fit, dashes: NLO one). 

In Tab. 14.41 we compare the values for as(Mz) obtained for the individual data sets at NLO and 
NNLO in the present fit and with results obtained by some of the experiments. Here the Qr.5-value 
for BCDMS was re-evaluated using the value of Ag^g = 224 MeV and a,(10 GeV) = 0.160 [l83l . 

These values correspond to a NLO fit with nf = 4 in the MS -scheme. We evolved this value back 
to the charm threshold keeping tif = 4 and determined then asiMz) evolving forward passing the 

bottom threshold. In [1591 higher twist contributions and A'^^T^ were fitted together for the BCDMS 

n f=4 

fip and fid data resulting in the somewhat larger value A^^^ = 263 + 42 MeV and the NLO value 
(^s(Mz) = 0.113 + 0.003 (exp), which was also obtained [iSOllllOL Both values are compatible 
within errors. 

In Fig. |4.18l we plot the -profile using eq. (lAlll) at NLO and NNLO. To that end, we compare 
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the fit result with the Qr^-behavior of the individual data set, fixing all other parameters. The 
minimum and variation (A;^'^ = 1) then determine the values given in Tab. 14.41 see also Fig. 14.181 
For BCDMS and NMC we find complete consistency to the values given by the experiments and 
the present analysis. The downward shift Aa^(Mz) which is consistently observed when going 
from NLO to NNLO amounts to values between 0.0030 and 0.0055, with a lower sensitivity for 
the DY data, which yield rather low values with large errors. The fitted central values are well 
covered by the individual data sets. Clearly, the above comparison is done ignoring the correlation 
to other parameters, and is therefore of more indicative character. The stability of the results, on 
the other hand, allows to conclude that fully compatible sets of precision data were used. 

We have also performed NLO fits, including the Tevatron jet data [I1051 - I108I . Furthermore, 
we have formally extended the analysis fitting the DIS and DY data at NNLO while treating the 
Tevatron jet data at NLO and supplementing threshold corrections based on soft gluon resum- 
mation HI 0911 for the single jet inclusive data. This latter approximation we denoted by NNLO*, 
cf. [I102L A NLO measurement of as(Mz) was also performed by CDF [I1L2J, with larger errors 
than in MM, a^iMz) = 0.1178 !°S^(exp) +2;j|;ji(scale) +0.0059 (PDF). At NLO the diff"erent 
sets of Tevatron jet data do not modify the value obtained in our standard analysis. A consistent 
NNLO is not yet possible since the corresponding scattering cross sections still have to be calcu- 
lated. Again a systematic downward shift of Aas{Mz) = 0.0029 - 0.0047 is obtained upon going 
from NLO to NNLO*. The corresponding central values are Icr compatible with our NNLO central 
value in eq. (14.131) . We would like to mention that already our former ABKM09 results [fT6l give 
a very good description of the CMS jet data [I113II and also the Tevatron 3-jet data [I114L We note 
that in a recent NLO analysis of the 5-jet cross section at LEP a value of as(Mz) = 0.1 156 ^^00034 
was obtained flll5L 




Figure 4.19: The correlation coefficient of a^iMz) with the nucleon twist-4 coefficients H2 (solid line) and 
Hj (dashes) versus x as obtained in our NNLO fit. 

The higher twist terms play an important role in the determination of from the DIS data ||59l . 
In our analysis they contribute up to 10% of the cross section at the low margin of and W given 
by eq. (13.31) . As a result, the value of Os is strongly correlated with the twist-4 coefficients, which 
are extracted from the fit simultaneously with as, cf. Fig. 14. 191 The correlation is more pronounced 
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Figure 4.20: The same as in Fig. l4.4l for the data points rejected in the analysis of II116I . 



for H2 at large x and for Ht at small x. The latter affects the determination of even in the case 
of a more stringent cut on W since the low-x part of the data is not sensitive to this cut. For the 
variant of our NNLO fit with the cut of eq. (|4.11l) imposed and the higher twist terms set to zero, 
we obtain the value of as(Mz) = 0.1191(6), much bigger than our nominal result of eq. (14.131) . 
For comparison, the same fit with the higher twist terms fixed at the values of Tab. 14.31 gives 
(Xs{Mz) = 0.1131(5) comparable with eq. (14.131) . Note that in both cases the error in as is smaller 
than one of eq. (|4.13l) despite the reduced data set used in the fit. This says, that the uncertainty in 
as is essentially controlled by the higher twist term variation. To get rid of the impact of the higher 
twist terms on as an even more stringent cut on is necessary, in addition to the cut on W. With 
the NNLO variant of our fit and using 



W^>12.5 GeV^. 



lOGeV^ 



(4.14) 



the value of as{Mz) = 0.1 134(8) is obtained, if the higher twist terms are set to zero and as{Mz) = 
0.1135(8), if the higher twist terms are fixed at the values of Tab. 14.31 From this comparisons we 
conclude that a? is pushed to larger values due to the neglect of the higher twist terms in the case 
of a cut as in eq. (|4.1 II) . which is commonly imposed in the global PDF fits. Likewise, it is less 
sensitive to the details of the fit ansatz. As we have found earlier Il27]| . the value of or 5 extracted 
from our fit is quite sensitive to the treatment of the NMC data Il88l . Normally, we use the NMC 
data on the cross section in the fit, cf. Sec. 13. II If, however, we employ instead the NMC data on F2 
extracted by the NMC collaboration with their own assumptions about the value of R, the value of 
as{Mz) increases by +0.0035 for the case of our earlier ABKM09 NNLO fit |l27l. In comparison, 
for the variant of the ABKM09 with the cut of eq. (14.1 II) imposed and the high-twist terms set to 
zero, the value of as{Mz) is shifted by +0.0003 only. This is in agreement with the size of the 
or, variation obtained in the variant of the MSTW08 fit with an improved treatment of the NMC 
data f28|. Note, however, that the details of the DIS data treatment employed in the improved 
analysis of [28], are still different from ours. In places, this has an impact on the as value. E.g., 
if we combine the errors in the NMC and HERA data in quadrature, as it is done in Il22ll28l . the 
NNLO ABKM09 value of as{Mz) is shifted upwards by +0.0029 and its sensitivity to the NMC 
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data treatment is reduced to +0.001 1 . This effect may in particular explain the relatively big value 
of asiMz) = 0. 1 164 observed in [|28ll with the cuts similar to ones of eq. (14.141) . 

The impact of the NMC data treatment on the fit has also been studied for NN21 in [I117II . 
where little effect on the PDFs has been found. However the value of R obtained in that analysis 
is still much lower than the one obtained in our fit and closer to the value of R used by NMC to 
extract the value of F2 from the cross section (compare Fig. 2 in [I117II and Fig. 1 in [|27l ). Note 
in this context that the value of x^/N DP « 1.7 obtained in [11171 for the NMC data is much bigger 
than in our case, cf.Tab. 14. 2[ Furthermore, in the NN21 study the value of a, is fixed. Therefore, 
its correlation with the PDFs is not considered. Finally, the expressions for cross sections used in 
the NN21 analysis do not include the power corrections in Q^, which are numerically important at 
small (compare eqs. (2), (3) in II117II with eq. (1) in (Ml, and eqs. (I3l1) . (IX2l) of the present 
article). These differences make a detailed comparison of the results of [27] and [[117[[ difficult. 

The value of a, obtained in our fit is substantially constrained by the BCDMS data of P82[[83L 
which almost entirely survive after the cut of eq. ([3.3[) . Meanwhile the authors of [,116[ suggested 
to cut in addition the most inaccurate BCDMS data with low inelasticity y. The value of as(Mz) 
reported in [[116[[ with such cut is by ~ 0.009 larger than the one obtained from the analysis of 
the whole set of the BCDMS data at NNLO. The pulls of the \ow-y data rejected in the analysis 
of [[116 1 with respect to our NNLO fit are given in Fig. [4.20[ The fit is in reasonable agreement with 
data within the errors. Furthermore, rejecting these data points from the NNLO fit, we obtain a 
value of as(Mz) = 0.1 139(12), which is somewhat bigger than the one in eq. ([4.13[) . The statistical 
significance of the shift, however, is marginal. The discrepant findings of [[116[[ concerning the 
impact of the low-y BCDMS data may appear due to the fact that the systematic uncertainties in 
the data are not taken into account in [1 16J. In our case the systematic errors are included into the 
value ofx^, cf. Sec. |Al Therefore the \ow-y data points with an enhanced systematic uncertainty 
have reduced weight and do not affect the fit. 



Experiment 


«.v(Mz) 


NLO,.,p 


NLO 


NNLO 


N^LO* 


BCDMS 

NMC 

SLAC 


0.1111+0.0018 
117 +0.011 

-0.016 


0.1138 + 0.0007 
0.1166 + 0.0039 
0.1147 + 0.0029 


0.1126 + 0.0007 
0.1153 + 0.0039 
0.1158 + 0.0033 


0.1128 + 0.0006 
0.1153 + 0.0035 
0.1152 + 0.0027 


BBG 
BB 




0.1148 + 0.0019 
0.1147 + 0.0021 


0.1134 + 0.0020 
0.1132 + 0.0022 


0.1141+0.0021 
0.1137 + 0.0022 



Table 4.6: Comparison of the values of a,(Mz) obtained by BCDMS [831 and NMC |[l03l at NLO with 
the results of the flavor non-singlet fits BBG [60] and BB [110] of the DIS flavor non-singlet world data, 
at NLO, NNLO, and N^LO* with the response of the individual data sets, combined for the experiments 
BCDMS [HaimiM, NMC [HI, and SLAC flM . 

We turn now to comparisons with other NNLO analyses which will be performed studying the 
contribution of different data sets to as{Mz). We first compare to the flavor non-singlet analy- 
ses [[60[[110[[ . The results are summarized in Tab. [4.6[ The values of as{Mz) at NLO turn out to be 
lower than those obtained in singlet analyses, cf. Tabs. [4.4| and [4.10[ However they are consistent 
within the scale variation errors. At NNLO both analyses lead to the same values. Also note the 
anti-correlation of the size of higher twist contributions in the large-x region with the inclusion 
of higher orders at leading twist, cf. [[60l [6n[67l[110L The next order, denoted by N^LO*, yields 
information on the remaining theoretical uncertainty. At N^LO*, the non- singlet three-loop Wilson 
coefficients are used [|9]JT0l and the four-loop non-singlet anomalous dimension is estimated with a 
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Pade-approximation and accounting for a 100 % error. In fact the latter extrapolation agrees within 
20% with the second moment of the non-singlet four-loop anomalous dimension II121[|122I . For 
the three experiments, which give the bulk information on as(Mz), the shift due to the N^LO* con- 
tributions amount to \Aas(Mz)\ = 0.0002 - 0.0006 and globally to 0.0007. At NNLO, the Qf,(Mz) 
values of the individual data sets vary by 0.0032, consistent within the Icr errors. 

Next, in Tab. 14.71 we compare with the fit results of NN21 ||123[I124II for individual data sets for 
DIS and other hadronic hard scattering data. At NLO the as{Mz) values range from 0.1 135 (E866, 
DY) to 0. 1 252 (NuTeV), with a corresponding range at NNLO of 0. 1 1 1 1 (DO jet) to 0. 1 225 (CDF 
jet). The value of a,(Mz) = 0.1204 + 0.0015 for the BCDMS data at NLO differs significantly from 
that given by the experiment 0.1111 +0.0018 [83]. Comparing the change of the a^(Mz) values 
between the NLO and NNLO analyses one finds downward shifts between 0.0075 (NuTeV) and 
0.003 (CDF R2KT) and upward shifts between 0.0055 (CDF Zrap) and 0.0061 (ZEUS H2), see 
Tab. 14.71 The values for as(Mz) obtained for the SLAC data are found to be larger than 0.124 
both at NLO and NNLO, cf. Fig. 2 in [124J. Thex^ values for the scans in ^^(Mz) at NNLO turn 
out to be worse than in NLO in the global analysis for the data sets of NMC, BCDMS, HERA L 
CHORUS, ZEUS F2C, DY E866, CDF Zrap, and DO Zrap, again see Fig. 2 in ITTM . Comparing 
the DIS only fit to the global analysis it is found that the values improve significantly, except 
for NMCp, and to a lesser extent for SLAC at higher values of q'.v(Mz), cf. Fig. 5 in [I124II . 



Experiment 


as(Mz) 




NLO,,p 


NLO 


NNLO 


BCDMS [82,831 


0.1111 


+ 0.0018 


0.1204 + 0.0015 


0.1158 + 0.0015 


NMCp [|84l 






0.1192 + 0.0018 


0.1150 + 0.0020 


NMCprf [|125| 
SLAC [|120| 


0.117 


+ 0.011 
- 0.016 


> 0.124 


0.1146 + 0.0107 
> 0.124 


HERA I [HI 






0.1223 + 0.0018 


0.1199 + 0.0019 


ZEUS H2 II1261I127H 






0.1170 + 0.0027 


0.1231+0.0030 


ZEUS F2C lil28H131ll 






0.1144 + 0.0060 




NuTeV |[95ll96l 






0.1252 + 0.0068 


0.1177 + 0.0039 


E605 [90] 






0.1168 + 0.0100 




E866 i91 [19211 13211 






0.1135 + 0.0029 




CDF Wasy il33[[ 






0.1181+0.0060 




CDF Zrap [[134|[ 






0.1150 + 0.0034 


0.1205 + 0.0081 


DO Zrap [jBSJl 






0.1227 + 0.0067 




CDF R2KT [105J 






0.1228 + 0.0021 


0.1225 + 0.0021 


DO R2C0N [fl07i 


0.1161 


+ 0.0041 
- 0.0048 


0.1141+0.0031 


0.1111+0.0029 


NN21 




0.1191+0.0006 


0.1173 + 0.0007 



Table 4.7: Comparison of the values of a,(Mz) obtained by BCDMS |[83l, NMC |[T03ll . and DO |[T04l at 
NLO with the results of NN21 II1231I124II for the fits to DIS and other hard scattering data at NLO and 
NNLO and the corresponding response of the different data sets analysed. 

In Tab. [M] the a, values determined by MSTW at NLO and NNLO [[I58II are compared for 
individual data sets in the fit. Here, the Icr errors, as defined by MSTW, are read off the corre- 
sponding plots in [1591. This definition of ' Icr' is obtained for values of Ax^ much larger than one. 
Moreover, these values do even strongly vary between the different measurements used, which is 
unlike the case for the ABM and NN21 analyses. This procedure leads to an enlargement of er- 
rors, which, e.g., in the case of BCDMS translates into a NNLO value as(Mz) = 0.1085 + 0.0095 
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Experiment 


«,v(Mz) 




NLO,;,p 


NLO 


NNLO 


BCDMSiup,F2 [[821 


0.1111 


+ 0.0018 


- 


0.1085 


+ 0.0095 


BCDMS iud,F2 [[831 






0.1135+0.0155 


0.1117 


+ 0.0093 


NMCjup,F2 [[Ml 


0.117 


+ 0.011 
- 0.016 


0.1275 + 0.0105 


0.1217 


+ 0.0077 


NMCjuJ,F2 [[841 




0.1265 + 0.0115 


0.1215 


+ 0.0070 


NMC jun/yup 11251 






0.1280 


0.1160 




E665//p,F2 [136l[ 






0.1203 




— 


E665 //J,F2 [136J 






— 




- 


SLACep,F2 [1621(891 






0.1 180 ±0.0060 


0.1140 


+ 0.0060 


SLACed,F2 [I62l[83 






0.1270 + 0.0090 


0.1220 


+ 0.0060 


NMC,BCDMS,SLAC, |[82[I841I120[| 






0.1285+0.0115 


0.1200 


+ 0.0060 


E886/NuSea pp, DY [|132[[ 






- 


0.1132 


+ 0.0088 


E886/NuSea pd/pp, DY [[911 






0.1173 + 0.107 


0.1140 


+ 0.0110 


NuTeVyA^.F9 [l 13711 






0.1207 + 0.0067 


0.1170 


+ 0.0060 


CHORUS vN,F2 11381 






0.1230 + 0.0110 


0.1150 


+ 0.0090 


NuTeVvA^,JcF3 [[137[l 






0.1270 + 0.0090 


0.1225 


+ 0.0075 


CHORUS vN,xF3 il38[[ 






0.1215 + 0.0105 


0.1185 


+ 0.0075 


CCFR [|951|%1 






0.1190 




— 


NuTeV vA^ ^ ///uZ [[95l[96l 






0.1150 + 0.0170 






HI en 97-00. cr^c [|86[[139H141II 






0.1250 + 0.0070 


0.1205 


+ 0055 


ZEUS 95-00. crf'^ n87[[142j-[144[| 






1235 + 0065 


0.1210 


+ 0.0060 


HI ec 99-00 cr'^'^ [[140[| 






0.1285+0.0225 


0.1270 


+ 0.0200 


ZEUS ep 99-00, cr^^ [| 14511 






0.1125 + 0.0195 


0.1165 


+ 0.0095 


Hl/ZEUS eo.FS'^^™ [|1281I129[I146W149[| 








0.1165 


+ 0.0095 


HI 99-00 incl. iets |1501 


0.1168 


+ 0.0049 
- 0.0034 


0.1127 + 0.0093 






7FTT^ /^n on inrl ipt« HI SI 111 S9II 




+ 0.0048 
- 0.0040 


U. 1 i / J + U.UUj J 






DO II pp incl. jets |107[[ 


0.1161 


+ 0.0041 
- 0.0048 


0.1185 + 0.0055 


1 1 33 


+ 0063 


CDF II ;?^incl. jets [|105|| 






0.1205 +0.0045 


0.1165 


+ 0.0025 


DU 11 W — > /V asvm. II1j3II 












CDF llW^lv asym. [Il54l 












DO II Z rap. P135[[ 






0.1125 + 0.0100 


0.1136 


+ 0.0084 


CDF II Z rap. l[155[l 






0.1160 + 0.0070 


0.1157 


+ 0.0067 


MSTW 




1202 + "■"''i^ 

^•^^^^ - 0.0015 


0.1171 


+ 0.0014 



Table 4.8: Comparison of the values of a,(Mz) obtained by BCDMS [[13, NMC HQM, HERA-jet |[T56l 
\T57 \. and DO f\M\ at NLO with the results of the MSTW fits to DIS and other hard scattering data at NLO 
and NNLO and the corresponding response of the different data sets analysed, cf. Figs. 7a and 7b in II158II . 
Entries not given correspond to a^iMz) central values below 0.110 or above 0.130; in case no errors are 
assigned these are larger than the bounds provided in form of the plots in Iil58[|159l . 



rather than the experimental one of Aas(Mz) = +0.0018. The latter accuracy is reflected in other 
analyses, however, cf. AasiMz) = +0.0007 (BB), AusiMz) = +0.0015 (NN21), and Aa,(Mz) = 
+0.0013 (ABMll). Similar effects are present for various other data sets, as can be seen com- 
paring the values given in Tabs. [4.4[ - [4.8[ In this way, almost all individual ai(Mz)-errors, even 
at NNLO, are larger than the typical theory uncertainty of about +0.0050 at NLO. We stress 



36 



Data Set 


ABM 11 


BBG 


NN21 


MSTW 


BCDMS 


0.1048 + 0.0013 


0.1126 + 0.0007 


0.1158 + 0.0015 


0.1101+0.0094 


NMC 


0.1152 + 0.0007 


0.1153 + 0.0039 


0.1150 + 0.0020 


0.1216 + 0.0074 


SLAC 


0.1128 + 0.0003 


0.1158 + 0.0034 


> 0.124 


( 0.1140 ± 0.0060 ep 
\ 0.1220 ± 0.0060 ed 


HERA 


0.1126 + 0.0002 




( 0.1199 ± 0.0019 
1 0.1231 ± 0.0030 


0.1208 + 0.0058 


DY 


0.101 +0.025 






0.1136 + 0.0100 




0.1134 + 0.0011 


0.1134 + 0.0020 


0.1173+0.0007 


0.1171+0.0014 



Table 4.9: Comparison of the pulls in a,(Mz) per data set between the ABMl 1, BBG JgOl, NN21 lfT24l 
and MSTW 11581 analyses at NNLO. 





as{Mz) 




BBG 

BB 

GRS 

ABKM 

JR 

JR 

ABMll 
MSTW 
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0.1134 + 0.0011 
0.1171+0.0014 
0.1173 + 0.0007 
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valence analysis, NNLO ff60B 
valence analysis, NNLO ||110|| 
valence analysis, NNLO ||160|| 
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HQ: BSMN-approach ^ 
dynamical approach [|20l 
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[|158| 
[|124| 
[|161| 


Gehrmann et al. 
Abbate et al. 


0.1153 + 0.0017+0.0023 
0.1135 + 0.0011+0.0006 


e^e~ thrust [162] 
e+e" thrust LI 6311 


3 jet rate 
Z-decay 
T decay 
T decay 
T decay 


0.1175 + 0.0025 
0.1189 + 0.0026 
0.1212 + 0.0019 
0.1204 + 0.0016 
0.1169 + 0.0025 


Dissertori et al. 2009 [11641 
BCK 2008/12 (N^LO) [I165N166II 
BCK 2008 11651 
Pich2011 [11671 
Boito etal. 2011 [[168[[ 


lattice 
lattice 
lattice 


0.1205 + 0.0010 
0.1184 + 0.0006 
0.1200 + 0.0014 


PACS-CS 2009 (2+1 fl.) 1169 1 

HPQCD 2010 [.170.1 

ETM 2012 (2+1 + 1 fl.) IITU 


BBG 
BB 


A 1141 +0.0020 
u.ii'+i _ Q 0022 

0.1137 + 0.0022 


valence analysis, N^LO(*) [I6OI 
valence analysis, N-^LO(*) [1 lOJ 


world average 


0.1184 + 0.0007 
0.1183+0.0010 


[172J (2009) 
[11671 (2011) 



Table 4.10: Summary of recent NNLO QCD analyses of the DIS world data, supplemented by related 
measurements using other processes. 

that the present analysis (ABMll) correctly accounts for the experimental systematic errors of 
all data sets used, cf. Tabs. 14.41 and 14.51 and that an enlargement of errors has not been neces- 
sary. In view of this fact, it is somewhat surprising that the final error in the MSTW analysis, 
i.e. AasiMz) = +0.0014 (MSTW), fully agrees with those obtained by BB, ABMll, and NN21, 
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cf. Tab. |4.9[ Unlike the case of NN21, the as(Mz) values of MSTW become generally lower at 
NNLO if compared to NLO, with the exception of the DO run II Z-boson rapidity data set, where 
the NNLO value is slightly higher. 

It is also interesting to compare the as{Mz) values obtained in the NN21 II123[|124II . MSTW 
II158II and ABMll analyses with respect to the individual data sets used in those fits. At NLO 
NN21 obtains lower values for DO R2C0N than MSTW and ABM and a significantly higher value 
for DO ZRAP than MSTW, see Tabs, gj] and |Ml At NNLO, the individual ff,(Mz) value for the 
data set CDF Zrap moves upward with respect to the NLO value with a significantly larger error, 
while for MSTW the value remains the same as at NLO. The values of ABMl 1 given in Tab. 14.51 
for Tevatron jet data are rather close to those of MSTW, both at NLO an NNLO*. For the NuTeV 
data the NLO and NNLO a^iMz) values show a bigger diff"erence for NN21 than for MSTW, while 
the NNLO values are rather similar. 

In Tab. |4]9] we compare the as values of the ABMll, BBG, NN21 and MSTW analyses for 
those data sets which are commonly used at NNLO. An NLO comparison would still be subject to 
a scale error of ~ 0.0050, which is usually to large to differentiate between the various fit results. 
For the BCDMS data ABMll, BBG and MSTW obtain lower values, while NN21 differs, e.g., 
by +4cr (or ^as{Mz) = 0.0057) from MSTW and +7(t (or Aa^Mz) = 0.01 10) from ABMl 1. For 
the NMC data ABM, BBG, and NN21 do agree very well, while the value of MSTW shows an 
upward shift of Aas{Mz) = 0.0066 compared to NN21. The size of both these shifts is of the order 
of NLO scale uncertainty and should not be present at NNLO. For the SLAC ep data the values 
by ABMll, BBG and MSTW are consistent within errors. On the other hand, MSTW reports a 
much larger value for the SLAC ed data than obtained by ABMl 1 and BBG. NN21 obtains partial 
as(Mz) values > 0.124 both at NLO and NNLO both in their global and DIS only analyses. This 
is in contrast to the present results, to BB, and to MSTW for the ep data. In the non-singlet BBG 
analysis the influence of the HERA data is strongly reduced, since most of these data are located 
within the quark-sea region. The fit results of NN21 and MSTW lead to values of as(Mz) ~ 0.120 
while those of ABMll yield a much lower value of 0.1126. Note that the MSTW analysis does 
not yet include the HERA run I combined data set [18]. The results for the DY data are consistent 
between ABMl 1 and MSTW, although the sensitivity of these data to as{Mz) is comparably small. 
In summary, despite the fact that NN21 and MSTW obtain nearly the same global fit values for 
(Xs{Mz), the above discussion shows that quite a series of individual pulls are different. Both 
ABMl 1 and MSTW do not confirm the relatively large as{Mz) values for the data sets of BCDMS 
and CDF R2KT and also the rise, at NNLO, of the NN21 value for CDF Zrap calls for further 
clarification, cf . Tab. 14. 7[ 

Summarizing the comparison of the present results on as{Mz) with the analyses ll60[|110[[T23l 
I124[|158[|159ll we observe a good agreement with ||60llllOII and find differences in a series of 
data sets for [[T23l[T24l[T58irT59ll both at NLO and NNLO, with partly different deviations in case 
of [1123012411 and II158[|159L NN21 does not agree with the BCDMS resuk. We find lower values of 
ocs{Mz) both for the HERA and the SLAC data. We do not confirm part of the as{Mz) values found 
in [|123lll24lll58lll59l for the jet data. We would like to mention once more that in our analysis 
no rescaling of errors is performed, which varies for different data sets in the analysis [I1581I159I , 
but we have accounted for the systematic errors given by the experiments directly. The procedure 
of l|158[|159l naturally leads to re- weighting of the impact of different data sets on the value of 
as{Mz). 

Finally we would like to summarize the results of different determinations of as{Mz) at NNLO 
and N^LO (or N^LO*) in Tab. I4.10[ Some part of these results has been reported in [I167L Flavor 
non-singlet analyses of the DIS world data were performed in ||60llllO[[160L with an accuracy 
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of AQr^(Mz) ^ 2% at NNLO. The difference between the value at N^LO* and NNLO amounts 
to ~ 0.0007, which provides an estimate for the size of a remaining uncertainty. The ABKM09 
analysis (T6\ is a combined flavor singlet/non-singlet fit of the DIS world data, DY and di-muon 
data. Here a remaining diff"erence of Aas(Mz) = 0.0006 due to the treatment of the heavy-flavor 
corrections was observed. These uncertainties signal the typical theory errors remaining at the 
present level of description. The JR09 analysis obtained very similar results in combined flavor 
singlet/non- singlet fits [20| with a slightly larger value in the standard fit compared to the dynami- 
cal approach. The present analysis obtains the same values as those found in [i l^l20l[60ll 11 Olll 6011 
before, while a slightly larger value asiMz) = 0.1 147 + 0.0012 was reported in ABMIO based 
on incorporating the combined HERA run I data [18|. However, the improved treatment of the 
heavy-quark contributions finally led to the present value as{Mz) = 0.1134 + 0.0011. The inclu- 
sion of Tevatron jet data, cf. [|l02l , although only a NNLO* analysis, alters this values at most to 
cts{Mz) = 0. 1 149 + 0.0012 and the effect of the complete, yet unknown NNLO corrections, remain 
to be seen. Low values of UsiMz) have not only been reported in analyses which are predominantly 
based on DIS data, but also from those of thrust in e''"e"-annihilation in [I163II , cf. also [I162II . 

Larger central values of as(Mz) at NNLO (and similar in size) are reported by MSTW II158II 
and NN21 II124II . These fits include a much broader set of hadronic scattering data in the analysis 
and above we have pointed out various differences between these analyses, see Tabs. 14.71 and 14. 8[ 
These differences manifest themselves in various cases in rather different pulls for the value of 
cts{Mz)- Note also that the Icr confidence level is very differently defined in the MSTW analysis 
compared to ABMll, NN21 and JR09. The (preliminary) central value of as{Mz) reported by 
CTIO [116111 is similar to MSTW and NN21 at NNLO, although accompanied by a rather large 
uncertainty of Aas(Mz) = 0.0050, which makes the CTIO result compatible with the lower values 
obtained in [|T6ll20l[60llll0lll60L We have also discussed the reasons for the difference between 
the central values of ar,(Mz) in the NLO and NNLO analyses |[l6l|20l[60l[IIOl[l6l and [I1241I158H 
and further comparisons may be performed. An earlier claim that this is caused by the Tevatron 
jet data is not confirmed, cf. ^U^. Larger central values for a^iMz) than in [|16[|20ll60llll0lll60ll 
are obtained for the 3-jet rate in e'^e" annihilation HI 6411 at NNLO and for the Z-decay width at 
N-^LO P1651, see Tab. 14.101 The current asiMz) values at NNLO extracted from r-decays vary 
between 0.1212 and 0.1 169 Ill65[|167l[r68] . One lattice measurement ni70 l| yields the same value 
as the current world average l|167[|172| . Other recent lattice results are compatible with this value 
and more lattice studies are still underway aiming at improved systematics 117311 . 

4.3 Heavy-quark masses 

The precise value of the heavy-quark masses is an important parameter in the description of DIS 
charm-quark production. In our fit we use the heavy-quark masses in the MS -scheme and their 
implementation for heavy-quark DIS discussed in [16] which allows us to relate the values for iric 
and Mb directly to PDG results [98J, as done in eq. (14.101) . 

The current DIS data displays great sensitivity to the charm mass as we have demonstrated 
previously ll32lll75L Therefore, based on the pseudo-data input from eq. (14.101) we have released 
the uncertainty of the quark masses to obtain the following results 

mdmc) = 1.27 + 0.06 GeV, mb(mb) = 4.19 + 0.13 GeV, (4.15) 

which shows that the error on from the inclusive DIS data used in the fit is comparable to the 
one quoted by the PDG ll98l (although not comparable to the single most precise measurement 
listed therein). Interestingly, we observe in the covariance matrix in Tabs. lAl 1"I - 1A131 correlations 
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Figure 4.21: Comparison of the data from 111741 for the semi-inclusive structure function different 
values of the momentum transfer versus x with predictions of various PDF sets at NLO and NNLO 
in QCD, all taken in the FFNS with = 3 and with a running-mass of m,. = 111 GeV iW], The NNLO 
predictions for F^' use the ABM 11 PDFs (solid curves), the JR09 PDFs [20] (dots), and the MSTW PDFs 
of (23 (long dashes). The NLO calculations are based on the NN21 PDFs of 1991 (short dashes). 



of nic with Us and some parameters of the gluon and the strange PDFs. The precision of DIS 
data to the value of has previously also been exploited for the first direct determination of the 
running mass for charm quarks from hadronic processes with space-like kinematics as a variant of 
ABKM09 yielding values consistent with but systematically somewhat lower than the PDG world 
average HBllfTTSll . 

It is interesting to compare the results of the fit to the most recent HERA data from the HI 
collaboration HI 741 for the charm structure function in heavy-quark DIS extracted with the 
use of the HVQDIS code II176II . This is done in Fig. 14.211 for our 3-flavor running-mass NNLO 
predictions with nic of eq. (14.151) using the approximate NNLO QCD predictions of ||39] - I4T1| . The 
data are not used in our fit, however the agreement is quite good, as well as for the predictions based 
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Figure 4.22: The charm- (left) and the bottom-quark (right) PDFs obtained in the global fit: The dotted 
(red) lines denote the ±lcr band of relative uncertainties (in percent) and the solid (red) line indicates the 
central prediction resulting from the fit with the running masses of eq. (I4.15I ). For comparison the shaded 
(grey) area represents the results of ABKM09 [,16il . 

on the NNLO 3-flavor PDFs of JR09 (201 and MSTW ES. The predictions using the NLO 3-flavor 
NN21 PDFs of [|99ll undershoot the data. The differences may be related to the peculiarities of the 
so-called general mass VFNS modeling employed in those fits [|22ll24ll99l . Also wrong analyses 
of the combined HI and ZEUS data on F'^ exist [11771 . 

Finally, we want to mention that the issue of heavy-quark masses also has consequences for 
heavy-quark PDFs, because the uncertainty on heavy-quark PDFs is directly related to the accuracy 
of the numerical values for the quark masses rric or m^. The latter appear parametrically in the 
OMEs used to generate charm- and bottom-PDFs in schemes with nf = 4 or nf = 5 flavors. With 
the results of eq. (|4.15l) and the use of the MS -scheme this uncertainty in the charm and bottom 
PDFs can be significantly reduced. 

In Fig. 14.221 we display the PDFs generated in this way. We obtain a charm-PDF with compa- 
rable uncertainties to the one of |fT6l (which then has used the pole mass definition for nic), while 
the resulting uncertainty of the bottom- PDF is greatly reduced, see also [I175II . This improvement 
has a significant impact on LHC phenomenology, as it will, e.g., allow for precise predictions for 
the production of single-top-quarks and other processes sensitive to bottom-PDFs. 

4.4 Moments of the PDFs 

In Tab. 14.1 II we summarize different values of the second moment of the valence quark densities 
obtained in NNLO analyses at the scale = 4 GeV^. It is evident, that these moments are rather 
stable quantities for all PDF sets considered as they are mostly influenced by the data normaliza- 
tion. They are closely related to the moments which are being measured in lattice simulations. Of 
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Table 4.11: Comparison of the second moment of the valence quark distributions at NNLO and N^LO 
obtained in different analyses at -4 GeV^. 
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Figure 4.23: Comparison of lattice computations for the second moment of the non-singlet distribution as 
a function of the pion mass m,^ with the result of ABM 11 given in Tab. |4.1l| along with the uncertainties of 
the respective measurement. 



central importance is the quantity 

1 



(xViQ^)) = Jdxx[[uix,Q^) + u,(x,Q^-[d(x,Q^) + ds{x,Q^}, (4.16) 







where q = qv + qs with q = u,d. 

In Fig. 14.231 the result for eq. (14.161) obtained in the present analysis is compared with recent 
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lattice computations using varying numbers of flavors as a function of the pion mass m„ employed 
on the lattice. In detail, these are QCDSF (n/ = 2) ETSl, RBC/UKQCD (n/ = 2 + 1) ([TTll, LHPC 
(n/ = 2 + 1) ETMC (n/ = 2) [181J and ETMC (n/ = 2 + 1 + 1) [|182[|183L It is apparent from 

Fig. 14.231 that there are substantial diff"erences, even for low pion masses, between those lattice 
measurements and the experimental determinations of Tab. |4~Tn 



5 Benchmarks for cross sections 

In this section we quantify the impact of the new PDF set on predictions for benchmark cross 
sections at the Tevatron and the LHC. To that end, we confine ourselves to (mostly) inclusive cross 
sections which are known to NNLO in QCD either completely or in very good approximation, 
see J26l for previous work along these lines. NNLO accuracy is actually the first instance, where 
meaningful statements about the residual theoretical uncertainty are possible, since at NLO the 
latter which is conventionally determined from a variation of the renormalization and factorization 
scale is generally still too large, given the precision of present collider data. 

In detail, we consider the following set of inclusive observables: hadronic W- and Z-boson pro- 
duction [[Tn[T2l . the cross section for Higgs boson production in the dominant channels, ggF lfT2l 
[TMMeH . vector-boson fusion (VBF) with VBFNNLO lfT87l[T88H . and in Higgs-strahlung [[T89l . 
The cross section for top-quark pair production is approximately NNLO (based on threshold re- 
summation, see e.g., [I190II ) and is computed with HATHOR [|191| . We also consider the lepton (/-) 
charge asymmetry in hadronic VF*-boson production as a function of the rapidity [fT3] - [r5l[T92| . 
Throughout the entire section we focus on the QCD corrections only. That is to say, we neglect all 
electroweak radiative effects at NLO, which often amount to corrections of (9(few)% at the LHC 
and, therefore, need to be considered in precision predictions. 

The PDF uncertainties quoted here are calculated by summing over the npDp sets provided by 
the various groups, where npoF is the number of parameters used in the fit. Typically, we quote 
the symmetric error according to 



A(TPDF = I (^o-o-fe)^ (5.1) 

where crk is obtained by using the k-th PDF f^, which parametrizes the +lcr- variation of the k-th 
fit parameter after diagonalization of the correlation matrix. In some cases, e.g., for MSTW Il22l . 
asymmetric PDF errors are provided, in which case the variation in the ^-th fit parameter is given 
by a pair of PDFs fl^'~. The resulting asymmetric PDF error is then computed according to 

^'^PDF = I max(0,-ha-yt,+ -cro,-Ha-yt--o-o)2, (5.2) 

^o-pjjp = I ^ mm{0,-o-k,+ + (To,-o-k-+(To)^ . (5.3) 

yk=l,npDF 

For MSTW, we are using the set with 68% confidence level error estimates throughout. 

In a Monte Carlo approach like the one advocated by NN21 [i23il . the PDF uncertainty can be 
determined as the quadratic deviation from the central fit as in Eq. (|5.1I) . but with an additional 
factor 1/ yfnpDF- For reasons of efficiency and run-times, we are using the NN21 PDF with 100 
sets in our comparisons only, see also the discussion in Sec. l5.5[ 
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5.1 W- and Z-boson production 



We start by presenting results for W- and Z-boson production at the LHC at = 7 TeV. For the 
electroweak parameters, we follow 126] and choose the scheme based on the set (Gf,Myy,Mz). Ac- 
cording to IMl, we have Gf = 1.16637 x 10"^ GeV'^, Mw = 80.399 + 0.023 GeV, Mz = 91.1876 + 
0.0021 GeV and the corresponding widths YiW^) = 2.085 + 0.042 GeV and r(Z) = 2.4952 + 
0.0023 GeV. The weak mixing angle is then a dependent quantity, with 

f = I ^ = 0.2307 + 0.0005 , (5.4) 

and p = 1.01047 + 0.00015. Finally, the Cabibbo angle Be yields the value sin^^c = 0.051. 

At NNLO the theoretical uncertainty due to scale variation is small compared to the PDF 
error, see Tabs. 15.11 and 15^ The change in the predictions between ABKM09 and ABM 11 is 
small. NN21 and MSTW typically predict smaller cross sections with differences at the level 
of l-2cr, while the numbers of JR09 are significantly smaller, see also the detailed discussion 
in [|26l. Most importantly, there is the choice to consider in particular the VK-boson cross section 
in alternative schemes with nf = A ox nf = 5 flavors, since contributions of the initial bottom PDFs 
being proportional to the CKM matrix element Vtt are kinematically suppressed. 
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Table 5.1: The total cross sections for gauge boson production at the LHC ( = 7 TeV) for different PDF 
sets and to NNLO accuracy The errors shown are the scale uncertainty based on the shifts jj. = M^jzl^ and 
ji = 2Mw/z and, respectively, the Icr PDF uncertainty. Numbers are in pb. 
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Table 5.2: Same as Tab.Ofor the PDF sets with «/ = 4. 



Comparing the results in Tabs . 15 . 1 1 and 15 . 21 for the PDF sets with = 4 and = 5 we observe 
that the numbers for the «/ = 5 scheme are always larger, the differences being less than Icr in the 
PDF uncertainty, though. These differences, which become successively smaller at higher orders, 
i.e. as we go from NLO to NNLO accuracy, originate from changes in the light flavor and the 
gluon PDFs when the bottom PDF is generated perturbatively, recall Sec. 14. 3[ In summary, the 
differences between the results in Tabs. 15.11 and 15.21 for a given PDF set constitute an intrinsic 
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uncertainty of the perturbative prediction. Comparisons of heavy-flavor PDFs including mass 
effects have also been studied in [I193I . 

Next, we address the charged-lepton asymmetry data I|1941I195| | as obtained by the ATLAS and 
CMS experiments and compare it to the NNLO predictions based on the ABMl 1 PDFs in Fig. l5.1[ 
All differential distributions for W- and Z-boson production are computed with the fully exclusive 
NNLO program DYNNLO [I151I192II . which allows to take into account the kinematical cuts imposed 
in the experiments (cf. Fig. 15 .11 ). see also [196] for an alternative code. 

The overall agreement with both experiments is sufficiently good, however at values of 77 ~ 1.5 
for the lepton pseudo-rapidity the data show a different trend with respect to the predictions. Pre- 
liminary data on the charge-lepton asymmetry at large rapidities obtained by the LHCb collabo- 
ration HI 9711 are also in good agreement with the ABKM09 predictions. To check the impact of 
the LHC charged-lepton asymmetry data on our fit we have performed a variant of the ABMll 
analysis which consists of adding the data of [ 1 94[|195L We have found, however, that the change 
in the PDF central values and their errors is only marginal in view of still big uncertainties in the 
data. 
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Figure 5.1: The data on charged-lepton asymmetry versus the lepton pseudo-rapidity r\ obtained by the AT- 
LAS [ 194 1 (left panel) and CMS II195II (right panel) experiments compared to the NNLO predictions based 
on the DYNNLO code 1115111921 and the ABMll NNLO PDFs with the shaded area showing the integration 
uncertainties. The ABKM09 NNLO predictions are given for comparison by dashes, without the integration 
uncertainties shown. 



5.2 Higgs boson production 

Let us now discuss the cross sections for the Standard Model Higgs boson production, where all 
dominant channels are known to NNLO in QCD. 

We start with ggF in Tabs. 15.3115.51 where the cross section is driven by the gluon luminosity 
and the value of a, from the effective vertex. The NNLO QCD corrections obtained in [fT2lll84| - 
II 8611 still lead to a sizable increase in the cross section at nominal values of the scale, i.e. /i = 



45 



m//. Further stabilization is achieved beyond NNLO on the basis of soft gluon resummation, see 
e.g., MM. 

We observe in Tabs. l5.3l - I5.5l that the ABM 11 predictions are rather stable with small changes 
only due to the gluon PDF discussed in Sec. IH The values of MSTW are typically larger than 
those of ABMl 1 and of ABKM09, roughly by (9(10%) depending on the Higgs mass and the LHC 
collision energy, which has direct consequences for the current Higgs searches at the LHC. E.g., at 
niH = 125 GeV MSTW predicts an 8% larger cross section of which 6.5% are due to the difference 
in as- In terms of PDF uncertainties, this discrepancy is significant at the level of 3-4cr. The 
reasons for the different gluon luminosities in the relevant .x:-range and the value of as(Mz) have 
been illustrated in Sec. |4l The great sensitivity of the ggF rate to constraints from higher orders in 
QCD in the treatment of fixed-target DIS data has already been discussed extensively in [|27ll . 
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ABM 11 



ABKM09 |fT6l| 



JR09 fionini 



MSTW [22J 



NN21 |[23 



100 
110 
115 
120 
125 
130 
140 
150 
160 
180 
200 
220 
260 
300 



21.31 
17.43 
15.85 
14.46 
13.23 
12.14 
10.30 
8.83 
7.63 
5.84 
4.58 
3.69 
2.55 
1.93 



+2.20 
-2.10 
+1.79 
-1.73 
+ 1.62 
-1.57 
+ 1.48 
-1.43 
+1.35 
-1.31 
+ 1.23 
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+ 1.04 
-1.03 
+0.89 
-0.88 
+0.76 
-0.77 
+0.58 
-0.59 
+0.45 
-0.46 
+0.36 
-0.37 
+0.25 
-0.26 
+0.19 
-0.20 



+0.46 
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+0.34 
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+0.32 
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+0.30 
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+0.28 
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+0.24 
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+0.21 
-0.21 
+0.19 
-0.19 
+0.15 
-0.15 
+0.13 
-0.13 
+0.11 
-0.11 
+0.08 
-0.08 
+0.07 
-0.07 



21.19 
17.31 
15.73 
14.34 
13.12 
12.04 
10.21 
8.75 
7.57 
5.79 
4.55 
3.67 
2.55 
1.95 



+2.21 
-2.11 
+ 1.79 
-1.72 
+ 1.62 
-1.57 
+ 1.47 
-1.43 
+ 1.34 
-1.31 
+ 1.23 
-1.20 
+ 1.03 
-1.02 
+0.88 
-0.87 
+0.76 
-0.76 
+0.57 
-0.58 
+0.45 
-0.46 
+0.36 
-0.37 
+0.25 
-0.26 
+0.19 
-0.20 



+0.60 
-0.60 
+0.49 
-0.49 
+0.45 
-0.45 
+0.42 
-0.42 
+0.38 
-0.38 
+0.35 
-0.35 
+0.31 
-0.31 
+0.27 
-0.27 
+0.24 
-0.24 
+0.19 
-0.19 
+0.16 
-0.16 
+0.14 
-0.14 
+0.10 
-0.10 
+0.09 
-0.09 



20.47 
16.91 
15.45 
14.18 
13.02 
12.01 
10.28 
8.90 
7.76 
6.04 
4.83 
3.97 
2.85 
2.23 



+ 1.99 
-1.87 
+ 1.63 
-1.54 
+ 1.48 
-1.40 
+ 1.35 
-1.28 
+ 1.24 
-1.17 
+ 1.14 
-1.07 
+0.97 
-0.92 
+0.83 
-0.79 
+0.72 
-0.69 
+0.56 
-0.53 
+0.44 
-0.42 
+0.36 
-0.35 
+0.25 
-0.25 
+0.19 
-0.20 



+0.70 
-0.70 
+0.56 
-0.56 
+0.50 
-0.50 
+0.45 
-0.45 
+0.41 
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+0.37 
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17.17 
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9.71 
8.44 
6.51 
5.17 
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2.97 
2.30 



+2.50 
-2.34 
+2.03 
-1.92 
+ 1.85 
-1.75 
+ 1.68 
-1.60 
+ 1.54 
-1.47 
+ 1.40 
-1.35 
+ 1.19 
-1.16 
+ 1.02 
-0.99 
+0.88 
-0.86 
+0.67 
-0.67 
+0.53 
-0.53 
+0.43 
-0.43 
+0.30 
-0.31 
+0.23 
-0.24 



+0.25 
-0.35 
+0.20 
-0.30 
+0.20 
-0.26 
+0.17 
-0.25 
+0.17 
-0.22 
+0.15 
-0.21 
+0.14 
-0.18 
+0.13 
-0.16 
+0.11 
-0.14 
+0.10 
-0.12 
+0.09 
-0.10 
+0.07 
-0.08 
+0.06 
-0.07 
+0.05 
-0.06 



24.16 
19.83 
18.07 
16.51 
15.14 
13.92 
11.86 
10.21 
8.86 
6.83 
5.42 
4.41 
3.11 
2.41 



+2.73 
-2.49 
+2.22 
-2.02 
+2.01 
-1.84 
+ 1.84 
-1.67 
+1.68 
-1.53 
+ 1.54 
-1.40 
+ 1.30 
-1.19 
+ 1.11 
-1.02 
+0.96 
-0.88 
+0.73 
-0.67 
+0.57 
-0.54 
+0.45 
-0.44 
+0.31 
-0.31 
+0.23 
-0.24 



+0.31 
-0.31 
+0.27 
-0.27 
+0.25 
-0.25 
+0.23 
-0.23 
+0.21 
-0.21 
+0.20 
-0.20 
+0.18 
-0.18 
+0.16 
-0.16 
+0.14 
-0.14 
+0.12 
-0.12 
+0.10 
-0.10 
+0.08 
-0.08 
+0.07 
-0.07 
+0.06 
-0.06 



Table 5.3: The total cross sections for Higgs production in ggF at the LHC {^/s = 1 TeV) for different PDF 
sets and to NNLO accuracy. The errors shown are the scale uncertainty based on the shifts /u = mn /2 and 
ji - 2mH and, respectively, the Icr PDF uncertainty. Numbers are in pb. 



Next in size comes the VBF channel. All numbers in Tab. 15. 61 are computed with the VBFNNLO 
program [|187[|188l in the structure function approach, which describes VBF as a double DIS 
process, where two (virtual) vector-bosons V, (independently) emitted from the hadronic initial 
states fuse into a Higgs boson. Although the structure function approach to VBF is not truly exact at 
NNLO it includes the bulk of the radiative corrections so that the remaining contributions, are both, 
parametrically small and kinematically suppressed, see 1188 1. The residual theory uncertainty 
based on the scale variation is rather small and the cross sections for all PDF sets considered in 
Tab. |5.6| agree well, typically within l-2cr of the PDF uncertainties. 

Last in line we consider the Higgs-strahlung process, that is the associated WH and ZH pro- 
duction using the NNLO QCD corrections of [189]. See also [199| for fully exclusive QCD cal- 
culations at NNLO. The dominant part of the hard partonic cross section is the same as for W- 
and Z-boson production discussed above in Sec. 15. 1[ so that essentially the same PDFs are probed. 
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rriH 


ABM 11 


ABKM09 jm 


JR09 1201I2T1 


MSTWC 


m 


NN21 |23| 


100 


26.91 


+2.72 
-2.57 


+0.56 
-0.56 


26.82 


+2.73 
-2.59 


+0.75 
-0.75 


25.64 


+2.44 
-2.28 


+0.91 
-0.91 


28.86 


+3.08 
-2.85 


+0.31 
-0.44 


30.35 


+3.36 
-3.03 


+0.37 
-0.37 


110 


22.17 


+2.21 
-2.13 


+0.46 
-0.46 


22.05 


+2.22 
-2.13 


+0.61 
-0.61 


21.31 


+2.01 
-1.88 


+0.72 
-0.72 


23.84 


+2.51 
-2.36 


+0.25 
-0.37 


25.07 


+2.74 
-2.48 


+0.32 
-0.32 


115 


20.22 


+2.02 
-1.94 


+0.42 
-0.42 


20.10 


+2.02 
-1.94 


+0.56 
-0.56 


19.52 


+ 1.83 
-1.71 


+0.65 
-0.65 


21.78 


+2.29 
-2.15 


+0.24 
-0.33 


22.91 


+2.49 
-2.26 


+0.29 
-0.29 


120 


18.51 


+ 1.84 
-1.78 


+0.39 
-0.39 


18.39 


+ 1.84 
-1.78 


+0.51 
-0.51 


17.96 


+ 1.68 
-1.58 


+0.59 
-0.59 


19.97 


+2.09 
-1.98 


+0.21 
-0.31 


21.00 


+2.28 
-2.06 


+0.28 
-0.28 


125 


16.99 


+ 1.69 
-1.63 


+0.37 
-0.37 


16.87 


+ 1.68 
-1.63 


+0.47 
-0.47 


16.53 


+ 1.54 
-1.44 


+0.53 
-0.53 


18.36 


+ 1.92 
-1.82 


+0.21 
-0.28 


19.30 


+2.09 
-1.89 


+0.26 
-0.26 


130 


15.64 


+1.55 
-1.51 


+0.34 
-0.34 


15.52 


+ 1.54 
-1.50 


+0.44 
-0.44 


15.29 


+ 1.42 
-1.33 


+0.48 
-0.48 


16.94 


+ 1.76 
-1.68 


+0.18 
-0.27 


17.80 


+ 1.92 
-1.74 


+0.24 
-0.24 


140 


13.36 


+ 1.31 
-1.29 


+0.30 
-0.30 


13.25 


+ 1.31 
-1.29 


+0.38 
-0.38 


13.16 


+ 1.21 
-1.14 


+0.41 
-0.41 


14.52 


+ 1.49 
-1.44 


+0.16 
-0.23 


15.25 


+ 1.63 
-1.48 


+0.21 
-0.21 


150 


11.51 


+ 1.13 
-1.12 


+0.26 
-0.26 


11.42 


+ 1.12 
-1.11 


+0.34 
-0.34 


11.45 


+ 1.05 
-0.99 


+0.35 
-0.35 


12.56 


+ 1.29 
-1.25 


+0.15 
-0.20 


13.19 


+1.40 
-1.28 


+0.19 
-0.19 


160 


10.01 


+0.98 
-0.97 


+0.23 
-0.23 


9.92 


+0.97 
-0.97 


+0.29 
-0.29 


10.03 


+0.91 
-0.86 


+0.30 
-0.30 


10.96 


+ 1.12 
-1.09 


+0.14 
-0.18 


11.51 


+ 1.22 
-1.10 


+0.17 
-0.17 


180 


7.74 


+0.75 
-0.76 


+0.19 
-0.19 


7.67 


+0.74 
-0.75 


+0.24 
-0.24 


7.88 


+0.71 
-0.67 


+0.24 
-0.24 


8.55 


+0.86 
-0.85 


+0.11 
-0.14 


8.96 


+0.94 
-0.85 


+0.14 
-0.14 


200 


6.15 


+0.59 


+0.16 


6.10 


+0.59 


+0.20 


6.36 


+0.56 


+0.20 


6.84 


+0.68 


+0.10 


7.18 


+0.73 


+0.12 


-0.60 


-0.16 


-0.60 


-0.20 


-0.54 


-0.20 


-0.68 


-0.12 


-0.69 


-0.12 


220 


5.00 


+0.48 
-0.49 


+0.14 
-0.14 


4.96 


+0.47 
-0.48 


+0.17 
-0.17 


5.26 


+0.46 
-0.45 


+0.17 
-0.17 


5.61 


+0.56 
-0.56 


+0.09 
-0.10 


5.89 


+0.59 
-0.57 


+0.10 
-0.10 


260 


3.53 


+0.34 
-0.35 


+0.11 
-0.11 


3.52 


+0.33 
-0.35 


+0.13 
-0.13 


3.84 


+0.32 
-0.33 


+0.14 
-0.14 


4.04 


+0.40 
-0.41 


+0.07 
-0.08 


4.23 


+0.41 
-0.41 


+0.08 
-0.08 


300 


2.72 


+0.26 
-0.27 


+0.09 
-0.09 


2.73 


+0.26 
-0.27 


+0.11 
-0.11 


3.06 


+0.25 
-0.27 


+0.13 
-0.13 


3.18 


+0.31 
-0.32 


+0.07 
-0.07 


3.33 


+0.31 
-0.33 


+0.07 
-0.07 



Table 5.4: Same as Tab.Ofor the LHC at = 8 TeV. 





ABM 11 


ABKM09 


lfT6l 


JR09 ||20l|2Tl 


MSTW |22| 


NN21 1231 


100 


66.79 


+6.13 
-5.63 


+ 1.31 
-1.31 


67.29 


+6.28 
-5.78 


+ 1.80 
-1.80 


62.23 


+5.46 
-4.92 


+2.62 
-2.62 


70.76 


+6.91 
-6.23 


+0.80 
-1.12 


74.18 


+7.50 
-6.54 


+0.78 
-0.78 


110 


56.35 


+5.11 
-4.77 


+ 1.06 
-1.06 


56.62 


+5.21 
-4.88 


+ 1.47 
-1.47 


52.77 


+4.58 
-4.14 


+2.11 
-2.11 


59.75 


+5.76 
-5.27 


+0.62 
-0.95 


62.63 


+6.26 
-5.47 


+0.66 
-0.66 


115 


52.01 


+4.70 
-4.40 


+0.99 
-0.99 


52.20 


+4.79 
-4.49 


+ 1.35 
-1.35 


48.82 


+4.21 
-3.80 


+ 1.92 
-1.92 


55.17 


+5.31 
-4.85 


+0.60 
-0.84 


57.86 


+5.74 
-5.04 


+0.62 
-0.62 


120 


48.14 


+4.35 
-4.07 


+0.93 
-0.93 


48.26 


+4.42 
-4.14 


+1.25 
-1.25 


45.33 


+3.89 
-3.53 


+ 1.75 
-1.75 


51.12 


+4.89 
-4.50 


+0.50 
-0.80 


53.58 


+5.31 
-4.64 


+0.58 
-0.58 


125 


44.68 


+4.02 
-3.78 


+0.85 
-0.85 


44.75 


+4.07 
-3.85 


+ 1.16 
-1.16 


42.13 


+3.60 
-3.26 


+ 1.59 
-1.59 


47.47 


+4.52 
-4.18 


+0.50 
-0.71 


49.77 


+4.91 
-4.30 


+0.54 
-0.54 


130 


41.59 


+3.72 
-3.53 


+0.80 
-0.80 


41.61 


+3.77 
-3.58 


+ 1.07 
-1.07 


39.32 


+3.35 
-3.02 


+ 1.45 
-1.45 


44.22 


+4.18 
-3.91 


+0.42 
-0.70 


46.35 


+4.55 
-3.99 


+0.51 
-0.51 


140 


36.28 


+3.23 
-3.09 


+0.70 
-0.70 


36.24 


+3.26 
-3.12 


+0.94 
-0.94 


34.46 


+2.90 
-2.65 


+ 1.22 
-1.22 


38.63 


+3.63 
-3.41 


+0.37 
-0.59 


40.49 


+3.94 
-3.47 


+0.45 
-0.45 


150 


31.92 


+2.82 
-2.72 


+0.61 
-0.61 


31.85 


+2.85 
-2.74 


+0.82 
-0.82 


30.49 


+2.56 
-2.32 


+ 1.04 
-1.04 


34.05 


+3.18 
-3.01 


+0.33 
-0.51 


35.69 


+3.45 
-3.04 


+0.41 
-0.41 


160 


28.32 


+2.49 
-2.42 


+0.54 
-0.54 


28.22 


+2.50 
-2.44 


+0.72 
-0.72 


27.16 


+2.25 
-2.06 


+0.90 
-0.90 


30.25 


+2.80 
-2.67 


+0.30 
-0.44 


31.69 


+3.04 
-2.67 


+0.37 
-0.37 


180 


22.75 


+ 1.98 
-1.96 


+0.44 
-0.44 


22.62 


+ 1.99 
-1.95 


+0.59 
-0.59 


22.03 


+ 1.80 
-1.65 


+0.69 
-0.69 


24.39 


+2.23 
-2.16 


+0.24 
-0.36 


25.56 


+2.41 
-2.15 


+0.31 
-0.31 


200 


18.74 


+ 1.62 
-1.61 


+0.37 
-0.37 


18.61 


+ 1.62 
-1.61 


+0.49 
-0.49 


18.33 


+ 1.47 
-1.38 


+0.54 
-0.54 


20.17 


+ 1.83 
-1.79 


+0.21 
-0.29 


21.16 


+ 1.94 
-1.79 


+0.27 
-0.27 


220 


15.78 


+ 1.36 
-1.36 


+0.32 
-0.32 


15.66 


+ 1.35 
-1.35 


+0.42 
-0.42 


15.61 


+ 1.23 
-1.17 


+0.45 
-0.45 


17.05 


+ 1.53 
-1.51 


+0.18 
-0.25 


17.90 


+ 1.61 
-1.52 


+0.24 
-0.24 


260 


11.91 


+ 1.01 
-1.04 


+0.26 
-0.26 


11.81 


+ 1.00 
-1.03 


+0.33 
-0.33 


12.03 


+0.91 
-0.91 


+0.34 
-0.34 


12.98 


+ 1.15 
-1.15 


+0.16 
-0.20 


13.64 


+ 1.18 
-1.16 


+0.19 
-0.19 


300 


9.80 


+0.83 
-0.86 


+0.23 
-0.23 


9.73 


+0.82 
-0.85 


+0.28 
-0.28 


10.11 


+0.75 
-0.77 


+0.29 
-0.29 


10.79 


+0.95 
-0.96 


+0.14 
-0.17 


11.34 


+0.95 
-0.97 


+0.17 
-0.17 



Table 5.5: Same as Tab.[53]for the LHC at - 14 TeV. 



although at slightly larger values of x. The numbers in Tabs. l5.7l and l5.8l do not contain the gluon 
induced contribution for o-(ZH) II200II . Top-quark mediated effects, which we neglect here, yield 
small perturbative corrections which are largely independent of the production model, that is to 
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ABM 11 



ABKM09 na 



JR09 1201I2T1 



MSTW 1 22 1 



NN21 lEa 



100 

no 

115 
120 
125 
130 
140 
150 
160 
180 
200 
220 
260 
300 



1.673 
1.513 
1.440 
1.373 
1.307 
1.244 
1.137 
1.038 
0.950 
0.800 
0.679 
0.580 
0.429 
0.324 



+0.024 
-0.022 
+0.025 
-0.019 
+0.021 
-0.021 
+0.022 
-0.018 
+0.023 
-0.018 
+0.025 
-0.014 
+0.016 
-0.015 
+0.019 
-0.013 
+0.018 
-0.011 
+0.016 
-0.009 
+0.013 
-0.008 
+0.011 
-0.006 
+0.010 
-0.003 
+0.008 
-0.002 



+0.020 
-0.020 
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-0.018 
+0.017 
-0.017 
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-0.016 
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-0.016 
+0.015 
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+0.014 
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+0.012 
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+0.010 
-0.010 
+0.009 
-0.009 
+0.008 
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+0.006 
-0.006 
+0.004 
-0.004 



1.643 
1.483 
1.411 
1.345 
1.285 
1.223 
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1.018 
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0.785 
0.669 
0.571 
0.424 
0.321 



+0.022 +0.012 
-0.026 -0.012 
+0.026 +0.011 
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-0.014 -0.011 
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-0.016 -0.010 
+0.020 +0.010 
-0.022 -0.010 
+0.022 +0.009 
-0.013 -0.009 



+0.021 
-0.013 
+0.021 
-0.011 
+0.017 
-0.012 
+0.016 
-0.007 
+0.012 
-0.007 
+0.012 
-0.005 
+0.009 
-0.003 
+0.007 
-0.003 



+0.008 
-0.008 
+0.008 
-0.008 
+0.007 
-0.007 
+0.006 
-0.006 
+0.005 
-0.005 
+0.005 
-0.005 
+0.004 
-0.004 
+0.003 
-0.003 



1.599 
1.453 
1.388 
1.321 
1.260 
1.203 
1.098 
1.006 
0.923 
0.781 
0.665 
0.570 
0.425 
0.323 



+0.023 
-0.017 
+0.014 
-0.020 
+0.016 
-0.022 
+0.018 
-0.018 
+0.017 
-0.020 
+0.017 
-0.017 
+0.018 
-0.016 
+0.015 
-0.014 
+0.016 
-0.014 
+0.013 
-0.010 
+0.011 
-0.008 
+0.010 
-0.006 
+0.008 
-0.004 
+0.006 
-0.003 



+0.022 
-0.022 
+0.021 
-0.021 
+0.020 
-0.020 
+0.019 
-0.019 
+0.019 
-0.019 
+0.018 
-0.018 
+0.017 
-0.017 
+0.016 
-0.016 
+0.015 
-0.015 
+0.014 
-0.014 
+0.012 
-0.012 
+0.011 
-0.011 
+0.009 
-0.009 
+0.007 
-0.007 



1.616 
1.460 
1.391 
1.324 
1.264 
1.203 
1.099 
1.003 
0.918 
0.773 
0.658 
0.562 
0.417 
0.316 



+0.025 
-0.034 
+0.025 
-0.027 
+0.026 
-0.028 
+0.025 
-0.023 
+0.023 
-0.023 
+0.026 
-0.021 
+0.024 
-0.019 
+0.021 
-0.016 
+0.018 
-0.013 
+0.019 
-0.011 
+0.014 
-0.010 
+0.013 
-0.007 
+0.010 
-0.005 
+0.008 
-0.003 



+0.029 
-0.029 
+0.026 
-0.026 
+0.025 
-0.025 
+0.024 
-0.024 
+0.023 
-0.023 
+0.022 
-0.022 
+0.020 
-0.020 
+0.019 
-0.019 
+0.017 
-0.017 
+0.015 
-0.015 
+0.013 
-0.013 
+0.011 
-0.011 
+0.008 
-0.008 
+0.007 
-0.007 



1.603 
1.448 
1.378 
1.318 
1.252 
1.195 
1.093 
0.999 
0.915 
0.771 
0.656 
0.561 
0.418 
0.316 



+0.028 +0.021 
-0.028 -0.021 
+0.027 +0.019 
-0.023 -0.019 
+0.029 +0.018 
-0.024 -0.018 



+0.021 
-0.028 
+0.026 
-0.019 
+0.027 
-0.018 
+0.021 
-0.020 
+0.020 
-0.019 
+0.017 
-0.014 
+0.018 
-0.011 
+0.016 
-0.009 
+0.012 
-0.007 
+0.010 
-0.005 
+0.007 
-0.003 



+0.017 
-0.017 
+0.016 
-0.016 
+0.016 
-0.016 
+0.015 
-0.015 
+0.013 
-0.013 
+0.012 
-0.012 
+0.011 
-0.011 
+0.009 
-0.009 
+0.008 
-0.008 
+0.006 
-0.006 
+0.005 
-0.005 



Table 5.6: The total VBF cross sections at the LHC {^fs = l TeV) for different PDF sets and to NNLO 
accuracy as computed with VBFNNLO l|187[|188l . Errors shown are the scale uncertainties evaluated by 
varying /i^ and/ij in the interval yU,->/^/ £ [2/4,42] and, respectively, the PDF uncertainties. Numbers are in 
pb. 
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Table 5.7: The total cross sections cr{WH) for associated Higgs production in the WH mode at the LHC 
{^/s = l TeV) for different PDF sets and to NNLO accuracy. The errors shown are the scale uncertainty based 
on the shifts fi - {mh + Mw)/2 and = 2(m// + M^) and, respectively, the Icr PDF uncertainty. Numbers 
are in fb. 



say of the parton luminosity, see also [2011 for a recent discussion. For the Higgs- strahlung the 
scale uncertainty is typically small and the predictions of the various PDF sets in Tabs. l5.7l and l5.8l 
agree well within the quoted PDF uncertainties. 
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Table 5.8: Same as Tab. 15 .7 1 for o-{ZH) and the scale uncertainty based on the shifts yu 
lu = 2{mH + Mz). 
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Table 5.9: The total cross section for top-quark pair production at the LHC {^/s = 7 TeV) for different PDF 
sets and to (approximate) NNLO accuracy as computed with HATHOR II191II as a function of the pole mass 
nit. The errors shown are the scale uncertainty based on the shifts iu = mjl and /u = Init and, respectively, 
the Icr PDF uncertainty. All rates are in pb. 
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Table 5.10: Same as Tab.[52]for the LHC at - 8 TeV. 



5.3 Top-quark pair production 

Here we present results for top-quark pair production at the LHC at ^/s = ^ TeV. The cross section 
is driven predominantly by the gluon luminosity and the value of as, much like in the case of Higgs 
production in ggF. By far the largest parametric dependence of the cross section resides, however, 
in the value of top-quark mass. This is currently quoted with an experimental uncertainty of less 
than 1% as = 172.9 +1.1 GeV based on the kinematic reconstruction from the decay products 
and comparison to Monte Carlo simulations [|98l . This value is commonly considered to be pole 
mass of the top-quark, although one should keep in mind, that the reconstruction of the top-quark 
momenta from the observed (color-neutral) hadron momenta carries a further intrinsic uncertainty 
of (9(Aqcd)- Top quark mass measurements in a well defined scheme for mt from the inclusive 
cross section at the Tevatron give lower masses, with a larger uncertainty though [|190[|202| . 

In Tab. |5.9| we summarize the cross section values for a range of top-quark masses 165 < < 
180 in the pole mass scheme. The theoretical uncertainty at (approximate) NNLO is quantified 
by a variation of scale in the range m,/2 < fi < 2m/. We see that the predictions of ABKM09 and 
ABMl 1 are largely the same, the latter being slightly smaller due to a smaller gluon PDF in the 
relevant x-range, cf. Fig. 14.121 In comparison to other PDFs, the predictions based on ABMl 1 in 
Tab. 15.91 are significantly smaller and it seems, that precision measurements of the cross section 
at the LHC can, potentially constrain the gluon PDF. This discriminating power, however, relies 
critically on the accurate knowledge of all other non-perturbative parameters, in particular m^. 

Cross section predictions using the running top-quark mass in the MS -scheme instead of the 
pole mass definition for in Tab. 15.91 generally display improved stability of the perturbative 
expansion and good properties of apparent convergence which is reflected in smaller uncertainties 
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nit 



ABM 11 



ABKM09 [161 



JR09 



mstwiEi 



NN21 |23| 



165 
166 
167 
168 
169 
170 
171 
172 
173 
174 
175 
176 
177 
178 
179 
180 



995.3 
969.5 
942.5 
916.3 
890.8 
866.4 
842.5 
819.5 
797.0 
775.7 
754.6 
734.5 
714.9 
696.0 
677.6 
659.7 



+ 10.9 

-43.8 

+8.7 

-44.8 

+8.1 

-43.5 

+7.9 

-42.4 

+7.8 

-41.1 

+7.6 

-40.4 

+7.4 

-39.0 

+7.1 

-37.9 

+7.3 

-37.0 

+6.8 

-36.2 

+7.0 

-34.8 

+6.7 

-34.3 

+6.6 

-33.3 

+6.4 

-32.2 

+6.3 

-31.5 

+6.2 

-30.8 



+28.4 
-28.4 
+27.7 
-27.7 
+27.0 
-27.0 
+26.2 
-26.2 
+25.7 
-25.7 
+25.0 
-25.0 
+24.4 
-24.4 
+23.8 
-23.8 
+23.2 
-23.2 
+22.8 
-22.8 
+22.1 
-22.1 
+21.7 
-21.7 
+21.2 
-21.2 
+20.7 
-20.7 
+20.2 
-20.2 
+ 19.7 
-19.7 



993.6 
962.7 
935.8 
910.0 
884.5 
860.5 
836.6 
814.1 
792.1 
770.8 
750.1 
729.9 
710.7 
691.8 
673.9 
656.2 



+ 10.5 

-43.3 

+ 13.1 

-39.0 

+ 12.8 

-37.9 

+ 12.5 

-37.0 

+ 12.6 

-35.8 

+ 12.2 

-34.9 

+ 12.0 

-33.6 

+ 11.5 

-33.1 

+ 11.3 

-32.1 

+ 11.1 

-31.1 

+ 10.9 

-30.2 

+ 10.9 

-29.2 

+ 10.2 

-28.8 

+ 10.2 

-28.0 

+9.6 

-27.7 

+9.6 

-26.8 



+35.4 
-35.4 
+34.5 
-34.5 
+33.6 
-33.6 
+32.9 
-32.9 
+32.0 
-32.0 
+31.2 
-31.2 
+30.4 
-30.4 
+29.6 
-29.6 
+28.9 
-28.9 
+28.3 
-28.3 
+27.6 
-27.6 
+27.1 
-27.1 
+26.4 
-26.4 
+25.8 
-25.8 
+25.1 
-25.1 
+24.5 
-24.5 



1076.9 
1043.1 
1014.8 
987.6 
960.8 
935.3 
910.0 
886.1 
862.7 
840.0 
817.9 
796.4 
775.8 
755.6 
736.9 
718.2 



+8.3 

-37.2 

+ 12.5 

-31.6 

+ 12.3 

-30.6 

+ 12.2 

-29.8 

+ 12.4 

-28.8 

+ 12.1 

-28.0 

+ 11.9 

-26.8 

+ 11.5 

-26.5 

+ 11.4 

-25.6 

+ 11.2 

-24.7 

+ 11.2 

-23.9 

+ 11.1 

-23.0 

+ 10.5 

-22.6 

+ 10.4 

-21.6 

+9.4 

-21.5 

+9.4 

-20.8 



+39.7 
-39.7 
+38.9 
-38.9 
+37.8 
-37.8 
+37.1 
-37.1 
+36.2 
-36.2 
+35.5 
-35.5 
+34.6 
-34.6 
+33.9 
-33.9 
+33.1 
-33.1 
+32.4 
-32.4 
+31.9 
-31.9 
+31.5 
-31.5 
+30.7 
-30.7 
+30.0 
-30.0 
+29.3 
-29.3 
+28.7 
-28.7 



1131.4 
1100.2 
1070.0 
1041.0 
1012.4 
985.3 
958.4 
933.1 
908.3 
884.2 
860.9 
838.2 
816.5 
795.1 
774.9 
754.9 



+ 11.6 

-50.5 

+ 11.1 

-49.1 

+ 10.8 

-47.7 

+ 10.6 

-46.6 

+ 10.8 

-45.1 

+ 10.5 

-44.0 

+ 10.4 

-42.5 

+9.9 

-41.8 

+9.8 

-40.5 

+9.6 

-39.3 

+9.5 

-38.2 

+9.5 

-37.0 

+8.8 

-36.5 

+8.8 

-35.4 

+8.2 

-35.0 

+8.3 

-33.9 



+ 18.4 
-20.3 
+ 18.0 
-19.9 
+ 17.6 
-19.3 
+ 17.1 
-18.9 
+ 16.7 
-18.4 
+ 16.3 
-17.9 
+ 15.9 
-17.5 
+ 15.6 
-17.1 
+ 15.2 
-16.7 
+ 14.9 
-16.3 
+ 14.6 
-16.0 
+ 14.3 
-15.6 
+ 14.0 
-15.2 
+ 13.7 
-14.9 
+ 13.3 
-14.5 
+ 13.0 
-14.2 



1174.6 
1141.9 
1110.3 
1080.0 
1050.5 
1022.0 
994.1 
967.7 
941.8 
916.7 
893.0 
869.9 
847.0 
824.9 
804.0 
783.1 



+ 14.0 

-46.0 

+ 13.4 

-44.3 

+ 13.7 

-42.8 

+ 13.1 

-41.6 

+ 13.0 

-40.5 

+ 12.9 

-39.4 

+ 12.8 

-37.8 

+ 12.3 

-37.0 

+ 11.8 

-35.8 

+ 11.8 

-34.6 

+ 11.2 

-34.1 

+ 10.4 

-33.8 

+ 10.2 

-32.8 

+ 10.0 

-31.8 

+9.2 

-31.4 

+9.3 

-30.6 



+20.8 
-20.8 
+20.3 
-20.3 
+ 19.8 
-19.8 
+ 19.3 
-19.3 
+ 19.0 
-19.0 
+ 18.5 
-18.5 
+ 18.0 
-18.0 
+ 17.6 
-17.6 
+ 17.1 
-17.1 
+ 16.6 
-16.6 
+ 16.3 
-16.3 
+ 16.0 
-16.0 
+ 15.6 
-15.6 
+ 15.2 
-15.2 
+ 14.9 
-14.9 
+ 14.6 
-14.6 



Table 5.11: Same as Tab.|l9]for the LHC at = 14 TeV. 



due the scale variation, see e.g. [I190II . The PDF uncertainties given in Tab. 15 .91 and the differences 
between the various PDF sets are, of course, largely unaffected. 

5.4 Hadronic jet production 

The hadronic jet production is sensitive to all nucleon PDFs, in particular to the large-x gluon 
distribution, which is not sufficiently sensitive to other processes employed in PDF fits. Therefore 
the inclusive jet production data obtained by the Tevatron collider experiments H 10511 lOVl are used 
in the PDF fits of [|22ll24ll99ll in order to provide better constraints on the large-x gluon. The 
calculation of the full NNLO QCD corrections to this process is still in progress (see 11203 [|204 ll 
and references therein). This precludes a consistent use of the Tevatron jet data in our NNLO PDF 
fit. Nevertheless in order to check any potential impact of the jet Tevatron data on our PDFs we 
have performed trial variants of the NNLO ABKM09 fit with the Tevatron jet data added ni02ll . 
The NLO QCD corrections [1205 [12061 and the partial (soft gluon enhanced) NNLO corrections due 
to threshold resummation [[109[[ have been computed with the FastNLO tool [[207[[208[[ . 

In order to allow more flexibility of the large-.x; gluon distribution we have added the term y2,g^^ 
to the polynomial of eq. ([4.9[) with an additional fitted parameter 72,g making sure that further 
expansions of this polynomial do not improve the quality of those fits. In general, the Tevatron jet 
data overshoot the ABKM09 predictions, nevertheless they can be smoothly accommodated in the 
fit. The typical value of x^/N DP ^ I is achieved with account of the error correlations for the jet 
data sets of [[105l - |108[[ once they are included into the NNLO ABKM09 fit. Meanwhile the various 
data sets demonstrate a somewhat different trend with respect to the ABKM09 predictions. E.g., 
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Figure 5.2: Cross section data for 1-jet inclusive production from the DO collaboration II107I as a function 
of the jet's transverse energy Et for the renormalization and factorization scales equal to Ej compared to the 
result of [161 (circles) and a re-fit including this data (squares) including the NNLO threshold resummation 
corrections to the jet production [.lOQJ . 



the off-set of the DO inclusive jet data [I107II does not depend on the jet energy Et and therefore 
may be attributed to the impact of the currently missing full NNLO corrections, cf. Figs. 15 .21 and 
I5.3[ In contrast, for the CDF data of II105II obtained with the kj jet algorithm the pulls rise with Et 
and can be reduced only by means of a modification of the PDF shapes. 



cr{H)[pb] 


ABKM09 


DO 1-jet inc. 


DO di-jet 


CDF 1-jet inc. 
(cone) 


CDF 1-jet inc. 

(kT) 


Tevatron(1.96) 
LHC(7) 


0.770(50) 
14.34(41) 


0.859(29) 
14.68(29) 


0.833(27) 
14.69(27) 


0.815(25) 
14.11(28) 


0.842(25) 
14.44(27) 



Table 5.12: The predicted cross sections for Higgs boson production in ggF with m// = 120 GeV at Tevatron 
( = 1.96 TeV) and at LHC ( ^ 7 TeV) from NNLO variants of the ABKM09 fit Ull corresponding to 
Tab. 14.51 The uncertainty in brackets refers to the Icr standard deviation for the combined uncertainty on the 
PDFs and the value of a^iMz). The values in bold correspond to the published result 1,26,1 . 



The values of as extracted from the trial ABKM09 fits with the Tevatron jet data included are 
compared with the nominal ABKM09 value in Tab. 14.51 At most, they are bigger by Icr, while for 
the CDF cone jet algorithm data [106J the central value of as is even the same. The predictions 
for the light Higgs production cross section, which are defined by the gluon distribution at jc < 0.1, 
are also not very sensitive to the constraints coming from the Tevatron data, cf. Tab. I5.12[ The 
impact of the Tevatron jet data on the large-x gluon distribution is more significant. However, in 
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Figure 5.3: Same as Fig. 15 .21 for the cross section data for 1-jet inclusive production from the CDF collab- 
oration using a kj jet algorithm II105I . 



this context we note that the trend of the first LHC data on the jet production with respect to the 
various PDF predictions is different from the Tevatron measurements. The ABKM09 predictions 
are in better agreement with the CMS and ATLAS inclusive data of HI 1 3112091 than the predictions 
based on the PDFs of [|22ll24ll99l . which were tuned to the Tevatron inclusive jet data. Jet data 
from LHC is still subject to large systematic errors, though. Note also that the Tevatron dijet and 
3-jet production data are in good agreement with the ABKM09 predictions [1210L in contrast to the 
case of inclusive jet production at Tevatron. These ambiguities in the data as well as the limitations 
in the current theoretical treatment prevent the use of hadronic jet data in our fit. 

5.5 LHAPDF library 

For the cross section computations presented here we have used the LHAPDF library [|211[|212| 
to interface to our PDFs and to those of other groups. For wider public use and to facilitate 
cross section computations by the interested reader, we provide the results of the current analysis 
in the form of data grids accessible with the most recent version lhapdf - 5.8.7 of the LHAPDF 
library, which can be obtained from http : / /projects .hepforge . org/lhapdf , Please follow 
the instructions in the LHAPDF package to install the LHAPDF library. 
In detail, we provide three NLO grids for = 3,4,5 flavors, 

abmll_3n_nlo.LHgrid (0+28), 
abinll_4n_nlo.LHgrid (0+28), 
abinll_5n_nlo.LHgrid (0+28), 
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with the central fit and 28 additional sets for the combined symmetric uncertainty on the PDFs and 
on Os from eqs. (|4.12l) . (14.131) . the heavy-quark masses from eq. (|4.15l) . and the deuteron correction, 
cf. Fig. l2.1[ Likewise at NNLO, we have 

abml l_3n_nnlo . LHgrid (0+28) , 
abml l_4n_nnlo . LHgrid (0+28) , 
abml l_5n_rmlo . LHgrid (0+28) . 

For future measurements of the strong coupling constant, e.g., from data on hadronic jet pro- 
duction as well as for detailed phenomenological studies of the parametric dependence of observ- 
ables on as, we also provide grids with fixed values of as, see the discussion in Sec. 14.21 At NLO 
these are 20 sets in the range a, = 0.1 1 .. .0.13 and at NNLO 16 sets for as = 0.105 . . .0.12, cf. also 
Fig. 14. 161 The sets for the as scan are denoted 

abml l_5n_as_nlo . LHgrid (Q+2Q) , 
abml l_5n_as_nnlo . LHgrid (0+16). 

During the computation of the cross sections with the LHAPDF library we have noticed signifi- 
cant differences in the run-time of our programs, when determining the cross section value and its 
uncertainty for a particular PDF set. Usually, this requires sampling of the error PDFs which de- 
scribe the error of the corresponding parameter in a given set. The run-time, of course, depends on 
the number of parameters in a given set, but also on the parametrization of the PDF grid and details 
of the look-up algorithm. Efficiency for phenomenological studies requires sufficiently precise and 
at same time fast cross section computations. 

Interestingly, for all observables discussed above, we have found the same pattern for the re- 
quired run-times. The grids of ABKM09 [16J and ABMll are both equally fast, i.e., ?abkmo9 - 
^^ABMU- Next come the JR09 [20l|21j grids which are slightly slower by roughly a factor ?jro9 - 
1.3 • ?ABMii followed by the grids of MSTW [|22J with ^mstw - 2 • ?abmii- The grids of NN21 [i21 
are last in this row with run-times up by roughly a factor of five, ?nn21 - 5 • ?abmii- This is cor- 
related with the enormous size of the latter grids of (9(100) MByte for NN21 with 100 PDF sets, 
as compared, e.g., to the size of 0(5) MByte for ABMl 1 or to (9(15) MByte for MSTW. Note that 
the NN21 grids with 1000 PDF sets have a size of about (9(1) GByte and typical run times up by 
another order of magnitude, making PDF error computations even more inefficient. 

6 Conclusions 

We have presented the PDF set ABMl 1 which determines the parton content of the nucleon and 
measures asiMz) at NNLO accuracy in QCD. In order to achieve a description of our analysis 
which is as complete and as transparent possible, we have provided an extensive discussion of all 
ingredients. This begins with the theoretical foundations which are most advanced and (almost) 
fully consistent at NNLO in QCD using well-defined renormalization schemes for all parameters 
involved. We have given detailed information on our treatment of the additional higher twist terms 
and nuclear corrections as well as on all data sets involved and on our fit ansatz. The results 
presented have been exposed to numerous consistency checks to test stability and the statistical 
quality. 

We have found good agreement with our previous results in ABKM09 and the PDF sets of 
ABMll are readily available with the LHAPDF library for precision phenomenology at the LHC 
and other hadron colliders. We have studied the differences with respect to other PDF groups using 
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the code OPENQCDRAD for standardized precision comparison. In this way some of the differences 
with respect to other PDF sets could be explained. The observed differences in the gluon PDFs are 
clearly outstanding in this context. This has important implications for the Higgs boson searches at 
the LHC and requires urgently further studies. In our present fit, thanks to OPENQCDRAD, we have 
provided a framework in which such questions can be addressed in the future. 

Despite the fact that various extractions of asiMz) have reached an impressive level of pre- 
cision of 1%, there exist still larger systematic differences which depend on both, the particular 
observable under consideration as well as on the specific analysis carried out for a given observ- 
ables (if compared to other analyses of the same observable). We have tried to analyse the status of 
the measurement of as{Mz) from world DIS data and the effect of including other hard scattering 
data. At present we feel, that not all data used in global analyses are of sufficient quality, consid- 
ering the precision currently envisaged for as(Mz). In the future, high quality data, possibly from 
the LHC, may improve this situation. 

It is therefore clear that the simplest way to reduce PDF uncertainties is not to discard non- 
global sets. Rather, it consists of very detailed studies of data sets, using the most complete 
theoretical descriptions, following D. Hilbert's request: "Wir mussen wissen - wir werden wis- 
sen!" [I213II . Precision QCD analyses will always refer to compatible sets of precise data with 
a detailed account of the systematic errors and all known theoretical corrections, which allows to 
measure all non-perturbative PDF parameters and the strong coupling constant as{Mz) as precisely 
as possible. 
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A Statistics 

A.l Statistical procedures 

In our analysis we infer the vector of fitted parameters 6 from the experimental data minimizing 
the;^f^ functional 

N 

x\0) = J]{fi{e)-ydE,j(fj(e)-yj), dAli) 

where fi{d) is the fitted model, yt are the measurements, and Ejj is the measurement error matrix 
with the indexes i,j running through all data points included into the fit. The error matrix Ejj 
is the inverse of the covariance matrix C,j. If the data are uncorrelated the co variance matrix is 
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diagonal and 

Qj = SijO-^a-j, (1A12) 

where cr,- are uncorrelated errors in the measurements y,. If, in addition, the data are subject to 
the point-to-point correlated systematic fluctuations, the off"-diagonal terms appear in the covari- 
ance matrix as well. For counting experiments the systematic errors are commonly multiplicative. 
Therefore the general form of the covariance matrix employed in our analysis reads 

Qj = mrijfi0)fj0) + 5ijCT,(Tj , (|A13) 

where is the relative correlated systematic error in the measurement ji. The uncorrelated errors 
may stem both from the statistical and systematic uncertainties. In the later case they are combined 
with the statistical ones in quadrature to obtain aj. The estimator based on the functional of 
eq. (|A[1I) is statistically efficient and asymptotically unbiased in the limit of A'^ ^ oo . In case of the 
analysis of correlated data the fitted parameters may be biased [I214L however for the definition 
of eq. (|A[3I) the estimator is nevertheless asymptotically unbiased [I215II . The errors in the fitted 
parameters given in Tab. 14. H and the correlations coefficients of Tabs. lAlTI - IA[3l are propagated from 
the uncertainties in the data with account of the correlations in the latter and correspond to the 
standard statistical criterion A;^^ = 1 . The errors in the data normalization factors for the selected 
experiments, cf. Tabs. |3.1| and [X2l which were fitted simultaneously with other parameters, are 
calculated in the same way. 

The detailed information on the experimental uncertainities in all data sets considered in the fit 
and listed in Sec. |3] are available from http : //arxiv . org as an attachment to the arXiv version 
of our paper. This includes in particular the systematic errors in the SLAC data [216J discussed 
in Sec. 13.11 and, likewise, the systematic errors in the NuTeV data and the corrections on the 
unmeasured phase space for the NuTeV and the CCFR data [I217II discussed in Sec. 13.31 

A.2 Parameter correlation matrix 

Here, we finally present the covariance matrix for the correlations of the fit parameters of ABMl 1 
discussed in Sec. HI cf. Tab. 14. H and eq. (14.101) for nic and m/,. The strong coupling as(Mz) is given 
in eq. (|4.13l) . while in Tabs. lA[l"l - IA[3l we quote the correlations for OsiMo) with /i^ = 1.5 GeV^. 
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